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The deepest hole that has ever been dug is about 12 km deep. Geochemists analyze samples from 
the Earth's crust and from the top of the mantle. Seismology can reconstruct the density profile 
throughout all Earth, but not its composition. In this respect, our planet is mainly unexplored. 
Geo-neutrinos, the antineutrinos from the progenies of U, Th and *'^K decays in the Earth, bring 
to the surface information from the whole planet, concerning its content of natural radioactive 
elements. Their detection can shed light on the sources of the terrestrial heat flow, on the present 
composition, and on the origins of the Earth. Geo-neutrinos represent a new probe of our planet, 
which can be exploited as a consequence of two fundamental advances that occurred in the last 
few years: the development of extremely low background neutrino detectors and the progress on 
understanding neutrino propagation. We review the status and the prospects of the field. 
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I. INTRODUCTION 

The deepest hole that has ever been dug is about 12 km deep, a mere dent in planetary terms. Geochemists 
analyze samples from the Earth's crust and from the top of the mantle. Seismology can reconstruct the density profile 
throughout all Earth, but not its composition. In this respect, our planet is mainly unexplored. 

Geo-neutrinos, antineutrinos from the progenies of U, Th, and K decays in the Earth, bring to Earth's surface 
information coming from the whole planet. Differently form other emissions of the planet (e.g., heat, noble gases), 
they are unique in that they can escape freely and instantaneously from Earth's interior. 

Detection of geo-neutrinos is becoming practical as a consequence of two fundamental advances that occurred in 
the last few years: a) development of extremely low background neutrino detectors a nd b) progress on understanding 
neutrino propagation. In fact, KamLAND has reported in 2005 ( Araki et al\ . l200"5ah evidence of a signal originating 



from geo-neutrinos, showing that the technique for geo-neutrino detection is now available. 

Geo-neutrinos look thus a promising new probe for the study of global properties of Earth and one has to examine 
their potential. Let us enumerate a few items which, at least in principle, can be addressed by means of geo-neutrinos^. 

What is the radiogenic contribution to terrestrial heat production? 

There are large uncertainties on Earth's energetics, both on the value of the heat flow (estimated between 30 
and 45 TW) and on the separate contributions to Earth's energy supply (radiogenic, gravitational, chemical 



^ Additional goals for geo-neutrinos {e.g., the distribution of radio-elements in the co re, discrimination among models of mantle circula- 
tion, and the possibility of detecting plumes in the mantle jFiorentini et aiX\20Q5<h ) appear presently too ambitious for the available 
technology. 
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. . . ). Estimates of radioactivity in the Earth's crust, based on observational data, account for at least some 
8 TW. The canonical Bulk Silicate Earth (BSE) model provides about 20 TW of radiogenic heat. However, 
on the grounds of available geochemical and/or geophysical data, one cannot exclude that radioactivity in the 
present Earth is enough to account for even the highest estimate of terrestrial heat flow. 

An unambiguous and observationally based determination of the radiogenic heat production would provide 
an important contribution for understanding Earth's energetics. It requires determining how much uranium, 
thorium and potassium are inside the Earth, quantities which are strictly related to the anti-neutrino luminosities 
from these elements. 

Test of the Bulk Silicate Earth model. 

The BSE model present s a chemical composition of the Earth similar to that of CI chondritic meteorites see. 



e.g., ( McDonoughl [2OO3I : iPalme and O'Neilll . [20031 ). The consistency between their composition and that of the 



solar photosphere points towards considering CI representatives of the material available in the pre-solar nebula 
and the basic material from which our planet has been formed. Some authors, however, have argued for a 
genetic relationship of our planet with other chondrites, such as enstatite chondrites, which are richer in long 
lived radioactive elements ( Javovl[l995t) . 



We remind that BSE is a basic geochemical paradigm consistent with most observational data, which however 
regard mostly the crust and an undetermined portion of the mantle. The global abundance of no element in 
the Earth can be estimated on the basis of observational data only. Geo-neutrinos could provide the first direct 
test of BSE (and/or its variants) by measuring the global abundances of natural heat radiogenic elements. 

Heat generating elements in the crust: a test of the estimated abundances. 

The amount of radioactivity in the Earth's crust is reasonably well constrained by observational data, with the 
exception of the lowest portion. Most of the uncertainty on the amount of radioactivity in the crust arises from 
the different estimates about the lower crust. In this respect, a detector located well in the middle of a continent, 
being most sensitive to geo-neutrinos from the crust, might provide a significant check of the estimates on the 
crustal content of heat generating elements. 

A measurement of heat generating elements in the mantle. 

The estimated content in the mantle is based on cosmochemical arguments and implies that abundances 
in de ep layers have to be much larger than those m easured in samples originating from the uppermost 
layer ( Jochum et all Il983t IZartman and HainesL [Ml). Uncertainties on the heat generating elements con- 



tent of the Earth essentially reflect the lack of observational data on the bulk of the mantle. A geo-neutrino 
detector located far from continents would be mainly sensitive to heat radiogenic elements in the whole mantle, 
as the oceanic crust is thin and poor in these elements. 

What can be said about the core? 

Geochemical arguments are against the presence of radioact ive elements in the c ore, although alternative hy- 
pothesis have been advanced see, e.g., (iHerndon. . il996 : RamaTlurthv et aLl , l2003l ). 

Present non directional detectors can say little about the core; however some extreme hypothesis can already be 
tested. If a natural fission reactor were present in the E arth's core, as advocated by Herndon in a series of paper 
( Herndonl . ll998l[2003l : lHerndon and IIollenbachll200ll ). it would produce antineutrinos with a spectrum similar 



to that of man-made reactors. An excess of "reactor like" antineutrinos events c ould be detected . A detailed 
analy sis already excludes a natural reactor producing more than about 20 TW ( Dve et all . 120061 : iFogli et~al\ . 
120051) . 

On the other hand, "non c'e rosa senza spine" ^. We list here the main difficulties and limitations encountered when 
detecting geo-neutrinos: 

• First of all, even huge detectors cannot provide more than some hundreds of geo-neutrino events per year. 

• Geo-neutrino events are to be disentangled from reactor neutrino events, which provide a severe background at 
many locations. 



^ There is no such thing as a rose without a thorn. 



4 



• Some 80% of the geo-neutrino events are expected to arise from uranium decay chain and only 20% from thorium 
chain. Due to the low yield, it will be hard to extract information on thorium abundance from the difference in 
the spectra. 

• Geo- neutrinos from K cannot be observed by means of inverse beta on free protons, the classical reaction for 
antineutrinos detection. 

• Present detectors cannot provide directional information. 

In the next section we shall outline the main properties (sources, spectra and cross sections) of geo-neutrinos and 
in section [Till we present how the field has evolved. Available information on the radioactivity content of the Earth 
is summarized in Sec. IIVI and the debated issue of the sources and flow of terrestrial heat is examined in Sec. [V] 
Section IVII presents a reference model for geo-neutrino production, i. e. a calculation of geo-neutrino fluxes based 
upon the best available information on Earth's interior. This model is refined in section IVIII for a specific location 
(the Kamioka mine, Japan) with a detailed calculation of the flux generated in the region. Section IVIIII provides a 
strategy for determining Earth's radioactivity from geo-neutrino measurements. This approach is developed in detail 
for KamLAND, the results of this experiment being presented and interpreted in section IIXI The role of reactor 
neutrinos, which are generally a significant background for geo-neutrino detection, is discussed in section [X] The 
prospects of the field are summarized in the final section. 

As a rule, when a section is divided into subsections, the first one contains an overview of the main points, so that 
the reader can decide whether the more detailed information presented in the foregoing subsections is of interest to 
him/her. 

II. GEO-NEUTRINO PROPERTIES 
A. Overview 

The natural radioactivity of present Earth arises mainly from the decay (chains) of nuclear isotopes with half-lives 
comparable to or longer than Earth's age^: ^ssy, '^^'^Th, *'^K, ^^^U, and ^'^Rb. 

Properties^ of these isotopes and of the (anti)neutrinos produced from their decay (chains) are summarized in 
Table HI Actually neutrinos are produced only in electron capture of ^°K. In contrast to the Sun, Earth shines 
essentially in antineutrinos. 

TABLE I Properties of ^^*U, ^^^Th, *°K, ^^^U, and ^''Rb and of their (anti)neutrinos. For each parent nucleus the table presents 
the natural isotopic mass abundance, half-life, antineutrino maximal energy (or neutrino energy), Q value, QeH = Q — (-5(1^,1^)), 
antineutrino and heat production rates for unit mass of the isotope (e^, sh), and for unit mass at natural isotopic composition 
(e^, e'h)- Note that antineutrinos with energy above threshold for inverse beta decay on free proton {Eth ~ 1.806 MeV) are 
produced only in the firsts two decay chains. 
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The energy of *^Rb neutrinos is so low that it is very unlikely that its flux could be measured. Also heat production 
from *^Rb is at the level of 1% of the total^. For these reasons, from now on we shall consider only U, Th, and 



^ Isotopes in the list have abundances and decay rates sufficiently large to give contributions of order 1% or more to the estimated 
radiogenic heat production: other radioactive elements such as ^'^^Lu, '^^'^Sm, ^*^Rn, give contribu tions of order 10~4 or less. 
In the Table and in the rest of the paper, unless differently specified, nuclear data are taken from ijFirestone and Shirlevl Il996h . 

^ This estimate is obtained assuming an abundance of **^Rb about 50 times the one of uranium. 
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and refer to these three elements as the Heat Generating Elements (HGEs) and to the antineutrinos from their decay 
(chains) as geo-neutrinos. 

For each isotope there is a strict connection between the geo-neutrino luminosity L (anti-neutrinos produced in the 
Earth per unit time), the radiogenic heat production rate Hr and the mass m of that isotope in the Earth: 

L = 7.46 X m(238u) + 31.94 x to(23^U) + 1.62 x m{^'^'^Th) + 23.16 x m(^°K) (1) 

Hr = 9.52 X m(238u) + 55.53 x m{^^^\]) + 2.67 x m{^^'^Th) + 2.85 x m(*°K) (2) 

where units are 10^4 s^\ 10^2 -y^ and 10" kg, respectively. By using the natural isotopic abundances in Table U these 
equations can be written in terms of the masses of the three elements^: 

L = 7.64 X m(U) + 1.62 x m(Th) + 27.10 x 10"^ x to(K) (3) 

= 9.85 X m(U) + 2.67 X m(Th) + 3.33 X 10"^ X m(K) . (4) 

The geo-neutrino spectrum depends on the shapes and rates of the individual decays, and on the abundances and 
spatial distribution of the terrestrial elements. It is shown in Fig. [T] for a specific model. 




0.5 1 1.5 2 2.5 3 3.5 

Energy [MeV] 



FIG. 1 Differential geo-neutrino luminosity, from (|Enomotd . |2005| ) °. One assumes the following global abundances: a('^^*U) = 
15 ppb, a(2='5U) = 0.1 ppb, ^(^^^Th) = 55 ppb, a(*"K) = 160 ppm l|McDonoughl. Il999l '). 

"Data are from Enomoto's web page: http://www.awa.tohoku.ac.jp/~sanshiro/geoneutrino/spectruiii/index.htnil 

The complete geo-neutrino spectrum depends on a large number of beta transitions in the uranium and thorium 
decay chains and it is essentially a result of theoretical calculations. These should be checked by measurements of the 
corresponding beta spectra, at least for the most important decays which contribute to the geo-neutrino signal: those 
of ^^"^Bi and ^■^'*Pam in the uranium chain, ^^^Bi and ^^^Ac in the thorium chain. 

Geo-neutrinos originating from different elements can be distinguished - at least in principle - due to their different 
energy spectra, e.5., geo-neutrinos with E > 2.25 MeV are produced only in the uranium chain. 

Geo-neutrinos from 238u g^j-^^j 232rp]_^ (^^^^ those from ^ssy g^j^j ^^K) are above threshold for the classical anti-neutrino 
detection reaction, the inverse beta on free protons: 

i>e+J5^e+-hn- 1.806 MeV . (5) 



^ The coefficients are slightly different from those quoted in (iFiorentini et al. I. l2005bl . l2003bh . which did not include 235 y contribution. 
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Note that anti-neutrinos from the Earth are not obscured by solar neutrinos, which cannot yield reaction ([5]). On 
the other hand , antincutrinos from nuclear power plants are a significant source of background, as first observed in 
( Laeag and discussed in more detail in section 1x1 

An order of magnitude estimate of the geo-neutrino luminosity can be obtained by assuming that a large fraction 
of the heat released from Earth, « 40 TW, arises from the decay chains of uranium and thorium. Table [T] shows 
that each of the N geo-neutrinos from each chain is associated with energy release AE « Q/N « 10 MeV, so that: 



L(V + Th) w H/AE w 2.5 x lO""^ s"' . (6) 

The order of magnitude of the produced flux is <i)(P™)(U + Th) « L/{At:R'^), where i?© is the Earth's radius. 
The flux arriving at detectors will be smaller than that produced due to neutrino oscillations, ^^^"''>(\J + Th) = 
(Pee)*^P™HU+ Th), where (Pee) « 0.6 is the average survival probability. All this gives: 

^(arr) (U + Th) W 2 X 10^ Cm^^ g-1 _ (7) 



This is a flux comparable to that of solar neutrinos from decay ( Castellani et a/.l . Il997^ ■ however the detection of 



geo-neutrinos is a much more difficult task: their smaller energy implies that the signal is smaller and is in an energy 
region where background is larger. 

For an order of magnitude estimate of the signal rate in a one-kton detector (containing some 10"^^ free protons), we 
observe that the cross section for inverse beta decay at few MeV is cr ~ 10"*'^ cm^ and the fraction of antineutrinos 
above threshold is / « 0.05. This gives a signal ^(U + Th) « cr/$('^") (U + Th)Np « 30 yr'^. 

More precisely, the signal rates S'(U) and 5(Th) in a detector containing Np free protons are 

'5(U) = 13 X — X — I- yr-i (8) 

^ ^ lO^cm-^s-i 1032 ^ ^ 

5(Th) = 4.0 X ^ l \ X yr-i , (9) 

^ ^ 106cm-2s-i 1032 ' ^ ^ 

where $("'^)(U) and $(^")(Th) are the fluxes of antincutrinos from 238^ g^j^^j ^pj^ arriving at the detector. 

Events rates are conveniently expressed in terms of a Terrestrial Neutrino Unit (TNU), defined as one event per 
10^2 target n uclei per year, or 3.17 x 10"'"' s~^ per target nucleus. This unit, which is analogous to the solar neutrino 
unit (SNU) (|Bahcal]lll989h . is practical since one kton of liquid scintillator contains about 10^2 free protons (the 
precise value depending on the chemical composition) and the exposure times are of order of a few years. 

Concerning the relative contributions of thorium and uranium to geo-neutrino events, Eqs. ([5]) and (O, together 
with Eq. ^ give: 

Sm $^-H^^^Th) ^ L(^32Th) ^ J_ ^(232Th) 

S{V) -"--^^^ $(arr)(238U) ^'^^ ^ ^238u) -le"" m(238U) " ^ 

Since one estimates that in our planet m(Th)/TO(U) « 4, one expects S'(Th)/S'(U) 1/4. Note that, although the 
global thorium mass is four times than that of uranium, it contributes just 1/5 of the total signal S{IJ + Th). 



B. Decay chains and geo-neutrino spectra from uranium and thorium 

One needs antineutrino spectra for two main reasons: the calculation of the specific elemental heat production and 
of the signal in the detector. 

Heat production rate is calculated by subtracting from the Q value the energy (E) of antineutrinos averaged over 
the whole spectrum. In the case of 238u g^^^ '^^'^Th chains^ the average antineutrino energy is about 8% and 6% of 
the total available energy: an error of 10% on the calculation of {E) is sufficient to determine the elemental heat 
production to better the 1%. For this reason in the literature individual determinations of beta spectra have not been 
used to determine neutrino energy loss. Instead, the approximate relationship that, on average, neutrinos carry 2/3 



''' Concerning ^'^^U, its contribution to the heat production is just a few per cent, so that the energy subtracted by antincutrinos is not 
relevant. 
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of the decay energy for beta decay has been applied ( van Schmus W.rI . Il995l ). This approximation can be checked 
or improved if the complete spectrum is known. 

For calculating the signal in a detector we need to integrate the spectrum times the cross section: only the spectrum 
above the detection threshold is needed for this aim. 

On these grounds we shall concentrate on the antineutrino energy spectra from ^ssy g^^^^ '^^'^Th decay chains. In 
general, the chain involves many different (3 decays and the total antineutrino spectrum results from the sum of the 
individual spectra. 

For each decay chain, if the sample of material contains rii nuclei of type i, the number of alpha and beta decays 
i ^ j per unit time is: 

Tij — riiXibij , (11) 

where A^, is the inverse of the mean-life and bi j is the branching ratio, bij ~ 1. The probability of each decay in 
the chain is: 

where h indicates the decay-chain head. The Ri,j form a network, with an isotope at each node. Generally the 
network has the following properties: 



— (by definition), 



• Rh.j = 1 (normalization); 

assuming that the chain is in secular equilibrium, one has: 

• X)fe Rk,i — J2j each node i (equilibrium). 

These three conditions fully determine the network^. 

In general the beta decay i j involves transitions to different nuclear states which yield spectra with different end- 
points: we call Iij;k the percentage intensity of the fc-th beta transition^ and fij-k^E) the corresponding antineutrino 
energy spectrum normalized to 1 (see below). 

Then the antineutrino spectrum generated from the sample is: 

f{E) = Y,R^.,J2^^'r,kkrME) ■ (13) 

ij k 

Lifetimes l/A^, branching ratios bij and intensities lijik, can be found in ( Firestone and Shir levl . [19961 ). 

A somehow delicate point is the expression to be used for the antineutrino spectra fij;k{E) of the /3 decay of nucleus 
i to the nucleus j into the state k. It can be derived from that for electron energy spectrum (f>ij-k{W) by using energy 
conservation: 

krME) = ^^,rMW)\w=w,^..-E ^ (14) 

where W is the total electron energy and Wmax = m^c^ -t- E^ax with Emax being the maximal neutrino energy for 
the transition and me the electron mass. 

For allowed decays the electron energy spectrum has the well-known universal shape: 

(b^,rAW)^^W{Wrnax-W)^W'-myr-'/^e^y\r{j + ty)\' , (15) 



where: 



W 



7=yT^(a^, y^aZ^==== , (16) 



It can be seen as a circuit where Rij are the currents and fe^j the inverse of the resistance, and where it flows a unit of current. 

Ou r notation corresponds to the normahzation 5I]fe^i,j;fc = 1- ^ different normaUzation, ^i,j;fc = is used 

in llFirestone and Shirlevl Il996l') . 
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with Z denoting the nuclear charge of the daughter nucleus and a the fine structure constant. iV is a normalization 
constant such that: 

dW4>^,rMi^) = ^ . (17) 

Equation p5|) is generally used to estimate geo-neutrino spectra and this requires a few comments. 

1) Equation (|15p considers the effect of the bare Coulomb field through the relativistic Fermi function. Electron 
screening and finite nuclear size effects are not considered. These provide corrections to the spectrum shape of 
order of few per cent, a quantity which is not significant in comparison with the uncertainties mentioned below. 

2) All important contributions actually arise from (first) parity forbidden decays. In this case the spectrum does not 
need to have a universal shape, since it involves also momentum-dependent nuclear matrix elements. However 
experimental data show that many forbidden decays of high-Z nuclei have spectra close to the allowed one: 
the theoretical explanation is that these decays are dominated by momentum-independent matrix elements or 
matrix elements whose relevant momentum is the electron momentum near the nucleus pi? k, Za, which is 
weakly dependent on the emerging momentum (^ approximation). This provides a partial justification for using 
equation (|15p . The resemblance with the allowed spectrum depends on the nucleus and it is difficult to study at 
low el ectron energy ; in few cases one finds significant differences^'', ^^''Bi (for an experimental review see, 
e.g., (|Danie]| . [19681) 1. 

3) Measurements of electron spectra would be very useful - in particular at low energy - in order to check the 
predictions for geo-neutrino spectra, which are mostly theoretical. In this respect an experimental study of the 
beta decay of ^^"^Bi would be most significant. 

Regarding the intensities hj-k, the experimental errors on some of them should be reduced: at the moment they 
imply a few percent uncertainty on the total geo-neutrino signal (see Tables IIIII and |V] and relative comments) . 



1. The ^^®U decay chain 

^•^^U decays into ^°^Pb through a chain of eight a decays and six (3 decays"'^"'^. In secular equilibrium the complete 
network (see Fig. [3]) includes nine /3-decaying nuclei^^ summarized in Table HIl 

TABLE II Beta decays in the ^''^U chain. For each decay we present the probability, the maximal antineutrino energy and 
the number of effective transitions, defined as those producing antineutrinos with E > 1806 keV. 
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Effective 








[kcV] 


transitions 


234 ^ 


^=«Pa 


1.0000 


199.08 





234p 

-Tarn 


_^ 234u 


0.9984 


2268.92 


1 


214pb ^ 


214Bi 


0.9998 


1024 





2"Bi ^ 


2"Po 


0.9998 


3272 


4 


2i0pb ^ 


210BJ 


1.0000 


63.5 





2iOBi ^ 


210pQ 


0.9999 


1162.1 





^s^Pa ^ 


.234u 


0.0016 


1247.15 





2i8po ^ 


21«At 


0.0002 


< 265 





206 rpj _^ 


206 pb 


0.0001 


1533.5 





210rpj _^ 


210pj3 


0.0002 


4391.3 


5 



Spectra of high-Z nuclei, that do not follow the allowed spectra, are explained theoretically by cancelations of dominant terms: a detailed 
knowledge of the relative weights and signs of the nuclear matrix elements becomes necessary. 

If we call Na the number of a decays and Nfj the number of /3 decays, A and Z {A' and Z') the atomic number and charge of the initial 
(final) nucleus, then Na = {A - A') /A and Nfj = Z' - Z + {A - A')/2 
This accounts for all branches with probability > 10~^. 
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FIG. 2 The ^38^ decay chain. The two nucUdes inside the grey boxes (^^*Pa and ^^''Bi) are the main sources of geo-neutrinos. 




FIG. 3 Goo-ncutrino spectra from the five main /3 decays of the U chain. All spectra are normalized to one decay of the 
head element of the chain. Since the ^•^*U chain contains six /3 decays, the integral from zero to the end point of the total 
spectrum is 6. Note that only 0.38 neutrinos per chain are above thresholds. 
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Only three nuclides (^'^'^Pa, ■^^^Bi, ^-'^^Tl) yield antineutrinos with energy larger than 1.806 MeV and contribute to 
the geo-neutrino signal. The contribution from ^lo^pj negligible, due to its small occurrence probability and the 
uranium contribution to the geo-neutrino signal comes from five (3 decays: one from ^"^^Pa and four from ^^"^Bi (see 
Table Hill and Fig. [3]). In fact, 98% of the uranium signal arises from the first two transitions in Tabic [nT] and an 
accuracy better than 1% is achieved by adding the third one. 

TABLE III Effective transitions in the ^^s-y gjjajjj jjj addition to quantities defined in Table ITll we present the intensity Ik, 
its error A/^, type and percentage contributions to the uranium geo-neutrino signal, and to the (U-|- Th) geo-neutrino signal. 
For this last column we assume the chondritic ratio for the masses (Th/U = 3.9), which implies that 79% of the geo-neutrino 
signal comes from uranium. 



i ~* j 


Ri,j 


[keV] 


Ik 




Type 
[%] 




S'u 
[%] 


Stot 


234pa„ ^ 234U 


0.9984 


2268.92 


0.9836 


0.002 


1st forbidden (0^ 


)-0+ 


39.62 


31.21 


214Bi ^ 214pQ 


0.9998 


3272.00 


0.182 


0.006 


1st forbidden 1~ 


^0+ 


58.21 


45.84 






2662.68 


0.017 


0.006 


1st forbidden 1~ 


^2+ 


1.98 


1.55 






1894.32 


0.0743 


0.0011 


1st forbidden 1~ 


^2+ 


0.18 


0.14 






1856.51 


0.0081 


0.0007 


1st forbidden 1^ 


^0+ 


0.01 


0.01 



In the last column of Table IIIII we show the contribution of each decay to the total (U -I- Th) geo-neutrino signal: 
this is calculated using a ratio of Th to U signal Sxh/'S'u — 0.270, that comes from the ratio between the average 
cross sections (cr)232xh/(o')23su = 0.127/0.404 = 0.314 (see subsection III.DI) and an assumed chondritic ratio^'^ for 
the masses TO(Th)/TO(U) = 3.9. Present errors on the intensities of the second and third decay of Table Hill imply 
corresponding errors to the total signal of 1.5% and 0.5%, respectively. 



2. The ^■^^Th decay chain 

232rpj^ decays into ^ospj-, through a chain of six a decays and four j3 decays. In secular equilibrium the complete 
network (see Fig. [4]) includes five /^-decaying nuclei^"' summarized in Table ITVl 

TABLE IV Beta decays in the 232rpj^ chain. For each decay we present the probability, the maximal antineutrino energy and 
the number of effective transitions, defined as those producing antineutrinos with E > 1806 keV. 



i j 


Ri,j 


Em, ax 


Effective 






[kcV] 


transitions 


228Ra ^ 228 


1.0000 


39.62 





22»AC ^ 


1.0000 


2069.24 


2 


212pb ^ 214Bi 


1.0000 


573.8 





2123. ^ 212p^ 


0.6406 


2254 


1 


208 rpj ^ 208 pj^ 


0.3594 


1803.26 






Only two nuclides (^^^Ac and ^^^Bi) yield antineutrinos with energy larger than 1.806 MeV. The thorium contri- 
bution to the geo-neutrino signal comes from three f3 decays: one from ^^^Bi and two from ^^^Ac (see Table Ivl and 
Fig. [5]). In fact, 99.8% of the signal arises from the first two transitions in Table fVl The present error on the intensity 
of the second decay of Table \V\ implies a corresponding error to the total signal of 0.9%. 



The corresponding ratio of fluxes is $Th/*u = (4/6) X (m(Th)/m(U)) X (238/232) X (ru/rTh) X (1/0.9927) = 0.8579. 
This accounts for aU the branches with probabihty > 10~^. 
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Half-life 




10.64 h I 



^ TIte Pb 



"°" 1d.299>is| 



FIG. 4 The ^^^Th decay chain. The two nuchdes inside the grey boxes (^^*Ac and ^^^Bi) are the main sources of geo-neutrinos. 




FIG. 5 Geo-neutrino spectra from the three main /3 decays of the ^^^Th chain. All spectra are normalized to one decay of the 
head clement of the chain. Since the ^-^^Th chain contains four (3 decays, the integral from zero to the end point of the total 
spectrum is 4. Note that only 0.15 neutrinos per chain are above thresholds. 
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TABLE V Effective transitions in the ^^-^Th chain. In addition to quantities defined in Table HVl we present the intensity Ik, its 
error A/^, type and percentage contributions to the thorium geo- neutrino signal, and to the total (U+ Th) geo- neutrino signal. 
For this last column we assume the chondritic ratio for the masses (Th/U — 3.9), which implies that 21% of the geo-neutrino 
signal comes from thorium. 



i - 


j 




Efnax 

[kcV] 


h 


A/fc 


Type 

[%] 




[%] 


•S'tot 


212Bi - 


^ 212po 


0.6406 


2254 


0.8658 


0.0016 


1st forbidden l'~' 


^0+ 


94.15 


20.00 


22»Ac - 




1.0000 


2069.24 


0.08 


0.06 


Allowed 3+ 


2+ 


5.66 


1.21 








1940.18 


0.008 


0.006 


Allowed 3+ 


4+ 


0.19 


0.04 



C. Geo-neutrinos from K 

""^K undergoes branching decay to ■^"Ca (via f3 decay) and ^°Ar (via electron capture), both of which are stable: 
the simphfied decay scheme of ''''K is shown in Fig. [6] The half life is 1.277 x 10^ yr, with a 10.7 percent probability 
of decaying to ''"Ar and an 89.3 percent probability of decaying to '*°Ca. All decays to ''"Ca proceed directly to the 
ground state, but most of the decays to ''''Ar reach an excited state, see Table IVll 



TABLE VI Decays of ^''K. For each decay we show the maximal antineutrino/neutrino energy, the intensity and the type of 
transition. 



I J 


[keV] 


Ik 


Type 








1311.09 


0.8928 


3rd forbidden 4^ 




0+ 




44.04 


0.1067 


1st forbidden 4~ 




2+ 


e- + «K ^ «Ar + ,^ 


1504.9 


0.00047 


3rd forbidden 4~ 




0+ 




482.9 


0.00001 


3rd forbidden 4~ 




0+ 



Kelley et al. ( Kellev et al\ . 19591) determined the sp ectrum of (3 particles emitted in the decay to *°Ca (Fig. [7]). 
From these data Van Schmus ( van Schmus W.R] , Il995l ) obtained a mean (3 energy of 0.598 MeV, or about 45% of the 
total; the remainder, 0.722 MeV (55%), is carried away by the antineutrino"'^^. 

We remind that the antineutrinos from [Emax — 1-311 MeV) are below the threshold for inverse beta on free 
protons. Note also that the monochromatic neutrinos from '^"K have a very small energy (44 keV). 



D. From cross sections to event rates 



As already mentioned, the classical process for detection of low energy antineutrinos is the inverse beta decay on 
free protons: 



Vf. + p e + n 



The threshold of the reaction is: 



thr _ [Mn+jnef^Ml _ 



2Af„ 



806 MeV 



The total cross section, neglecting terms of order Ee/Mp, is given by the standard formula: 

a = 0.0952 X ( x lO-^^^m^ 
MeV^ 



(18) 



(19) 



(20) 



We checked that by using Eq. UlSt times the non-relativistic correction factor appropriate for a 3rd forbidden decay, S(p^,p^) 



-Pg + Tpipiipf, +pi), one finds the same average energy as Van Schmus. Note that the value of the maximal energy used by Van 



Schmus, Wm ax 
0.588 MeV. 



1.32, should be replaced with the more resent value: Wn 



1.31109. In this case the average /3 energy becomes 
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FIG. 6 Simplified decay scheme for K. 




FIG. 7 Experimental spectrum of electron kinetic energy for the decay of ^''K into ^"Ca, from (jKellev et all 1 19591 ). The 
circles show the measured spectrum including background, 1.46 MeV gamma and finite resolution corrections. The x-s show 
the spectrum after the electron escape corrections. The flags represent total estimated error at each point, due mostly to the 
uncertainty in the electron escape correction. 



14 



where Eg = Eg^ — [Mn — Mp)c^ is the positron energy, when the (small) neutron recoil is neglected, and Pf. is the 
corresponding momentum. The numerical factor in Eq. ([^0]) is tied directly, see (Bemporad et ai, 2002), to the 
neutron lifetime, known to 0.1% ( Yao et o/l . l2006| l. This expression of the total cross section is shown in Fig. [H 




2.5 3 
Ev [MeV] 



FIG. 8 Total cross section for + p ^ + n as a function of the antineutrino energy, Eq. (f20 



Corrections to the cross section of order E^/Mp, which are negligible for geo- neutrinos whereas should be con- 
sider ed at reactor energ i es, an d the angular distribution of the po sitrons are described by Vogel and Beacom in 
Ref. ( Vogel and Beacoml . Il999( ): see also Ref. (jBemporad et aLl . l2002f ). 

A more general dis cussion of the neutrino /nuc leon cross-section for energies from threshold up to several hundred 
MeV can be found in f Stru mia and Vissanil . [20031 ): in the same paper Strumia and Vissani give a simple approximation 
which agrees with their full result within few per-mille for E^, < 300 MeV, 



aiDeP) ~ 10-«[cm^]pe£^eS; 



-0.07056+0.02018 In B^-0. 001953 In-^ E„ 



Ep — E,, 



(21) 



where A is the neutron-proton mass difference and all energies are in MeV. They conservatively estimate an overall 
uncertainty of the cross section at low energy of 0.4%. At the energy relevant for geo-neutrinos (< 3.27 MeV) Eq. ((20|) 
overestimates the full result of Strumia and Vissani by less than 1% and it is already identical at about 2 MeV. 
The geo-neutrino event rate from the decay chain of element X = ^ssy 232 rpj-^ jg. 



S{X) = Npj dE^e{Ei,)a[E^)c^'^^"\Ei,) 



(22) 



where Np is the number of free protons in the target, e is the detection efficiency, (t{Ep) is the cross section for 
reaction (HH]), and: 



dr 



p{r) ax{r)C_ 



X 



12 Txmx 



fx{E,)p (^E„\R-r\ 



(23) 



is the differential flux of antincutrinos from 238jj qj. 232rpj^ arriving into the detector, p is the density, ax is the 
elemental mass abundance, Cx, "Of, and mx are the isotopic concentration, lifetime and mass of nucleus X. The 
energy distribution of antincutrinos fx{Ep) is normalized to the number of antincutrinos nx emitted per decay chain: 



nx 



dE^fx{E^) 



(24) 



p (^E[,, \ R ^ rlj is the survival probability for v with energy Ep produced at r to reach the detector at R. 

In yiew of the values of th e oscillation length one can average the survival probability over a short distance, see 
Ref. ( Mantovani et al] . \2004) . and bring out of the integral the averaged survival probability: 
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In this way we are left with: 

S{X)^Np{Pee) I dE^e{E,)a{E,)fx{E,) I dr J'x{r)Cx ^^e) 

J Jvs, 47r|i?-r|2 Txmx 

The second integral is proportional to the (angle integrated) produced flux of anti-neutrinos 

HX)^p^f dfP^^^ . (27) 

Note that this quantity is different from the flux normal to earth surface. Note also that "produced" essentially means 
the flux which one would observe in the absence of oscillations. 

One can also assume the detection efficiency as approximately constant over the small (< 2 MeV) energy integration 
region. Then Eq. ([26]) becomes 



S{X) = N,{P,,)e^{X) [ d£;,^i:^£lM^ . (28) 
J nx 

It can be useful to introduce an average cross section: 

{<y)x = J dEMEu)fx{E^)/ J dE-,fx{E-,) . (29) 

This is computed by using equation (j20p for the cross section a{E[,) and the spectrum fxiEp) obtained in the previous 
section. Thus one finds (cr)238u = 0.404 x lO^'^'^cm^ and ((7)232^11 = 0.127 x lO^'^'^cm^. 
The event number can thus be written as the product of a few terms: 

Six) = Np{P,,)enX){a)x (30) 

The result is 

5(-U) = 4.04 X 10- s- X (^) ( JSZ'^.L ) (31) 

(32) 



V1032 J Vl06cm-2s-i 
5(-Th) = 1.27 X 10- s- X (P.)s (^) 



"cm ^s 



This is the way in which Eqs. ^ and ^ were derived. Our goal in the rest of the paper will be to provide calculations 
of the produced fluxes based on geological models. 

It is interesting to examine the differential geo-neutrino signal per unit flux as a function of the energy: 

^^=<y{E,)fx{E,)/ JdE,fx{E,) . (33) 

This quantity is shown in Figs. [9l and 1101 for uranium and thorium, respectively. Note that most of the geo-neutrino 
flux originates from very few transitions. 



III. A HISTORICAL PERSPECTIVE 

Geo-neutrinos have been conceived during the very flrst attempts of neutrino detection, performed at the Han- 
ford nuclear reactor by Reines and Cowan in 1953. Experimental results showed an unexpected and unexplained 
backgrounds^. While on board of the Santa Fe Chief Train, Georg Ganiow wrote to Fred Reines (see Fig. [TT|) : 

It just occurred to me that your background may just be coming from high energy beta-decaying 
members of U and Th families in the crust of the Earth. 

The first estimate of geo-neutrino flux was given in a teletype message by Reines (Fig. [T^ in response to the letter of 
Gamow: 



Actually the background was due to cosmic radiation. 



Bi [1 .86 MeV] 




1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

Energy [MeV] 



FIG. 9 Geo-neutrino diflerential signal per unit flux from the five main /3 decays of the ^-^^U chain, see Eq. (|33p. 
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"V-ri^)^* ti TKA o.^ W On lay** 



FIG. 11 The message from Georg Gamow to Fred Raines. 



Ml MMM It lot KHK OH as lunis, WIUVUtIM HUI 

ft 

M» tJKAT TBI Un ORLX DOOOt BIB»X tOB AMW U» , l| Hw 

■iinntrPtni ct^ tm sac. m» ib lor n io? oi m. 



FIG. 12 The teletype message from Reines to Gamow. 
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Heat loss from Earth's surface is 50 erg cm ^ s ^. If assume all due to beta decay than have only 
enough energy for about 10^ one-MeV neutrinos per cm^ and s. 



In the scientific literature, geo-neutrinos were introduced by Eder ( Eded . [T966| l in the sixties and Marx ( Marxl . 
Il969l ) soon realized their relevance. In the eightie s Krau ss et al. discussed their potential as probes of the Earth's 
interior in an extensive publication |Kra uss et al). '1984'). In the nineties the first paper on a geophysical journal 
was published by Koba yashi et al. (jKobayash i and F ukao, 1991). Of particular interest, in 1998, Raghavan et al. 
(jRaghavan et aLl . ll998t) and Rothschild et al. ()Rothschild et al\ . 
Borexino for geo-neutrino detection. 



1998t) pointed out the potential of KamLAND and 



In the last few ye ars more papers appeared than in the previous decades: in a series of papers ([Fiorentini et al\ . 
l2003al |2004 L2005b"'c'd], l2003bl : iMantovani et all l2004 f) Fiorentini et al. discussed the role of geo-neutrinos for de- 
termining the radiogenic contribution to the terrestrial heat flow and for discriminating among different models of 
Earth's composition and origin. A reference model for geo-neutrino production, b ased on a compos i tional map of 
the Earth's crust and on geochemical modeling of the mant le, was presented in (fM antovani et al*.. '2004'). bimi- 
lar calculat ions were per f ormed by Enomoto et al. (Enomo to et al\ . l2005i) and by Fogli et al. (^ogli_et _aZ., 2006,). 
The claim (jEguchi et of an in dication of geo-neut r ino events in the first data, release of KamLAND stim- 

ulated several theoretical investigat ions toomogatsk v et "al. . 200d lEguchi et al}. '2003'; 'Fields an d Hochmuthl. l2006t 
iFiorentini et all l2005al : iFogh all . 2005; McKeown and Vogell . 12004: Miramonti, 2003; Nunok awa et all. l2003l). A 
summary of the theoretical predictions is presented in Fig.[T31 Early models (|Edeii . ll966l : lKobavashi and Fukaol . ll99ll : 
iMarxl Il969; ) (full circles) assumed a uniform uranium distribution in the Earth and different values of the uranium 
mass. In fact these predictions are almost proportional to the estimated mass of heat generating elements. The huge 
signals predicted by Eder and by Marx were obtained by assuming that the uranium density in the whole Earth is 
about the same as that observed in the continental crust; Marx (Eder) assumed thus an uranium mass 30 (60) times 
larger than that estimated within the BSE model (see Section |V.C|) . 
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FIG. 13 Previous estimates of the geo-neutrino signal S, renormalized to the average survival probability (Pee) = 0.59, and the 
corres ponding estimated uran ium mass m. The signal is in Terrestrial Neutrino Units (1 TNU = 1 event /year/ 10"^^ proton). 
From ([Fiorentini et al ]. l2005dl ). 

Krauss et al. ([Krauss et al\ . 1 19841) distributed about 10^^ kg of uranium uniformly over a 30 km crust. The 
other estimates (crosses) are all obtained by using the BSE value for the uranium mass (« 10^^ kg) as an input 
and differ ent mod els for distributing the uranium content between crust and mantle. In this class, Rothschild et al. 
|Rothschild et al\ . |1998 ) obtained the minimal prediction by assuming for the crust a very small uranium abundance, 
definitely lower than the values reported in more recent and detailed estimates. 

In July 2005 the KamLAND collaboration presented the first evidence of a signal truly originati ng from geo- 



neutrinos, showing that the technology for geo-neutrino detection is now available. KamLAND reported ( Araki et al\ . 
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l2005al ) data from an exposure of Np = (0.346 ±0.017) x 10"^^ free protons over a time of 749 days. In the energy region 
where geo-neutrinos are expected, there are 152 counts. After subtracting several backgrounds, there remain about 25 
true geo-neutrino events. This indicates the difficulties of this experiment: a signal rate of one geo-neutrino event per 
month, to be distinguished over a five times larger background, mostly originating from the surroundi ng nuclear power 
plants . The implication of KamLAND result on radiogenic terrestrial heat have been discussed in ( Fiorentini et all . 
l2005af ). 

Following the important KamLAND result, a meeting speci fically devoted to stu dy the potential of geo-neutrinos 
in Earth's science was gathered at Hawaii in December 2005 ()Learned et aLl . [2"006t) . It provided a first opportunity 
for a joint discussion between the communities of particle physics and of geo-science. 

In a few years KamLAND should provide definite evidence of the geo-neutrino signal, after accumulating a much 
larger statistics and reducing background. In the meanwhile other projects for geo-neutrino detection are being 
developed. Borexino at Gra n Sasso, which is expect ed to take data soon, will benefit from the absence of nearby 
reactors. Domogatski et al. ( Domogatskv et aLl , |2006[ ) are proposing a one-kton scintillator detector in Baksan, again 
very far from nuclear reactors. A group at the Sudbury Neutrino Observatory in Canad a is studying the possibility 
of using liquid scintillator after the physics program with heavy water is completed ( Chenl . l2006l ). The LENA 
proposal envi sages a 30-kton Hquid scintillator detector at the Center for Underground Physics in the Pyhasalmi mine 
(Finland) (U ndagoitia et aLl . l2006t ). Due to the huge mass, it should collect se veral hundreds o f event s per year. The 
proposal of a geo-neutrino directional detector at Curacao has been advanced ( de Meiier et aLl . [20061 ). The possibihty 
of a detector located at Hawaii islands has been presented by Dye (|Dve et al\ . 20061 ). In conclusion, one can expect 
that within ten years the geo-neutrino signal from uranium and thorium will be measured at a few points on the 
globe. 



IV. RADIOACTIVITY IN THE EARTH 
A. A first look at Earth's Interior 

A global look at Earth's interior is useful before entering a detailed discussion on the element distributions. The 
amount of information which we (assume to) have on Earth's interior is somehow surprising, if one considers that the 
deepest hole which has ever been dug is only about twelve kilometers deep. 

Seismology has shown that Earth is divided into several layers, which can be distinguished from discontinuities in 
the sound speed, see Figs. [T3] and [TS] The outer layer is the relatively thin crust which accounts for 0.47% of the 
Earth mass; it is divided in two types, continental crust (CC) and oceanic crust (OC). The former averages 38 km in 
thickness, varying around the globe from 20 to 70 km, and it is made primarily of light elements such as potassium, 
sodium, silicon, calcium, and aluminum silicates. The oceanic crust is much thinner, from about 6 to 8 km. 

Inside this crustal skin is Earth's mantle which is 2900 km deep over all. Largely made up of iron and magnesium 
silicates, the mantle as a whole accounts for about 68% of Earth's mass. One distinguishes the upper mantle (UM) 
from the lower mantle (LM), however, the seismic discontinuities between the two parts do not necessarily divide the 
mantle into layers. The main questions about the mantle are: does it move as a single layer or as multiple layers? 
Is it homogeneous in composition or heterogeneous? How d oes it convect? The se questions sound simple, but the 
answers are complex, possibly leading to more questions, see ( Davies et aLl . l200"^ . 



Inside the mantle is Earth's core, which accounts for about 32% of Earth's mass. Based on comparison with 
the behavior of iron at high pressures and temperatures in laboratory experiments, on the seismic properties of the 
core, and on the fact that iron is the only sufficiently abundant heavy element in the universe, the core is generally 
believed to be made primarily of iron with small amounts of nickel and other elements. Over thirty years ago, 
however, it was suggested that a significant amount of potassium could be hidden in Earth's core, thus providing a 
lar ge fraction of the terrestria l heat flow through "^ ^K decay. Th i s cont roversial possibility has been revived recently 
in ( Rama Murthv et al\ . \2QQ'^ . see, however, Ref. ( Corgne et aLl . [20071 ). 



Concerning the density profi le of our planet, a classical r eferen ce is the Preliminary Reference Earth Model (PREM) 
of Dziewonski and Anderson ( Dziewonski and Anderson[ . Il98ll ). This one-dimensional spherically symmetric model 



is at the basis of all calculations for geo-neutrino production from the mantle. In the last twenty years seismic 
tomography has progressed so as to provide three dimensional views of the mantle. Density differences with respect 
to the one-dimensional model (typically of order of few percent) are most important for understanding mass circulation 
inside the mantle; however, they are too small in order to affect significantly the calculated geo-neutrino production. 

From seismic studies one can derive the density profile of our planet and the aggregation state of the different layers; 
however, one cannot reconstruct its composition. 

Earth global composition is generally estimated from that of CI chondritic meteorites by using geochcmical argu- 
ments which account for loss and fractionation during planet formation. Along these lines the Bulk Silicate Earth 
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FIG. 15 PREM (Preliminary Reference Earth Model) fPziew onski and AndersonL Il98lh velocity structure trough the Earth: 
p — density, a — seismic P-waves velocity, /3 — S-wavcs velocities " . 

"Figure taken from http: //shadow. eas . gatech . edu/'anewman/ classes/geodynamics/random/prem_earth.pdf 
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(BSE) model is built, which describes the element composition of the "primitive mantle", i.e. the outer portion of 
the Earth after core separation and before the differentiation between crust and mantle, see Table IVlIl The model is 
believed to describe the present crust plus mantle system. It provides the total amounts of U, Th, and K in the Ea rth, 
as these lithophile elements should be absent in the core. Estimates from different authors jM cDonoughl, 2003) are 
conco rdant within 10-15%, extensive reviews being provided in Refs. ( McDonough and Sunl .[T995: Palme a nd O'Neilll . 
I2003D . From the mass, the present radiogenic heat production rate and neutrino luminosity can be immediately 
calculated by means of Eqs. Q, (l2|), ([3]) and ([T]), and are shown in the Table IVlTTl 



TABLE VII The composition of the sil icate Earth. Abun dances are given in /ig ■ g ^ (ppm), unless stated as "%" which are 
given in weight percentage. Data from (jMcDonoughl . [200^ ) . 



H 


100 


Zn 


55 


Pr 


0.25 


Li 


1.6 


Ga 


4 


Nd 


1.25 


Be 


0.07 


Ge 


1.1 


Sm 


0.41 


B 


0.3 


As 


0.05 


Eu 


0.15 


C 


120 


Se 


0.075 


Gd 


0.54 


N 


2 


Br 


0.05 


Th 


0.1 


0(%) 


44 


Rb 


0.6 


Dy 


0.67 


F 


15 


Sr 


20 


Ho 


0.15 


Na (% ) 


0.27 


Y 


4.3 


Er 


0.44 


Mg (%) 


22.8 


Zr 


10.5 


Tm 


0.068 


Al (%) 


2.35 


Nb 


0.66 


Yb 


0.44 


Si (%) 


21 


Mo 


0.05 


Lu 


0.068 


P 


90 


Ru 


0.005 


Hf 


0.28 


S 


250 


Rh 


0.001 


Ta 


0.037 


CI 


17 


Pd 


0.004 


W 


0.029 


K 


240 


Ag 


0.008 


Re 


0.0003 


Ca (% ) 


2.53 


Cd 


0.04 


Os 


0.003 


Sc 


16 


In 


0.01 


Ir 


0.003 


Ti 


1200 


Sn 


0.13 


Pt 


0.007 


V 


82 


Sb 


0.006 


Au 


0.001 


Cr 


2625 


Te 


0.012 


Hg 


0.01 


Mn 


1045 


I 


0.01 


Tl 


0.004 


Fe (%) 


6.26 


Cs 


0.021 


Pb 


0.15 


Co 


105 


Ba 


6.6 


Bi 


0.003 


Ni 


1960 


La 


0.65 


Th 


0.08 


Cu 


30 


Ce 


1.68 


U 


0.02 



TABLE VIII Mass, heat production, and geo-neutrino luminosity of U, Th, K according to BSE. 





m 

[10^^ kg] 


He 
[10^2 2] 


Li, 
[10^* s-i] 


u 


0.8 


8.0 


6.2 


Th 


3.2 


8.7 


5.3 


«K 


1.1 


3.6 


26.5 



The BSE is a fundamental geochemical paradigm. It is consistent with most observations, which however regard the 
crust and the uppermost portion of the mantle only. Its prediction for the present radiogenic production is 19 TW. 

Concerning the distribution of heat generating elements, estimates for uranium in the (continental) crust based on 
observational data are in the range: 



mc{U) ^ (0.3^0.4) X 10" kg 



(34) 
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The crust - really a tiny envelope - should thus contain about one half of uranium in the Earth. For the mantle, 
observational data are scarce and restricted to the uppermost part, so the best estimate for its uranium content 
wm(C^) is obtained by subtracting the crust contribution from the BSE estimate: 

ruM {U) = niBSE [U)-mc[U) . (35) 

We remark that this estimate is essentially based on a cosmo-chemical argument and there is no direct observation 
capable of telling how much uranium is in the mantle, and thus on the whole Earth. 

Similar considerations hold for thorium and potassium, the relative mass abundance with respect to uranium being 
globally estimated as a(Th) : a(U) : a(K) « 4 : 1 : 12000. Geochemical arguments are against the presence 
of radioactive elements in the (co mpletely unexplored ) core, as discussed by McDonough in an excellent review of 
compositional models of the Earth (jMcDonoughl . |2003[ ). 

A comprehensive reyiew about the knowledge of Earth's interior is given in volumes 2 and 3 
of ([Holland and Turekianl . [20031 ). 



The following subsections are devoted to present, in some more detail, the available information on the amounts of 
heat generating elements in the whole Earth and within its separate layers. 



B. The BSE model and heat generating elements in the interior of the Earth 

In the BSE frame, the amount of heat/neutrino generating material inside Earth is determined through the following 
steps: 

(a) From the compositional study of selected samples emerging from the mantle, after correcting for the effects 
of partial melting, one establishes the absolute primitive abundances in major elements with refractory and 
lithophile character, i.e. elements with high condensation temperature (so that they do not escape in the 
processes leading to Earth formation) and which do not enter the metallic core. In this way primitive absolute 
abundances of elements such as Al, Ca and Ti are determined, a factor about 2.8 times CI chondritic abundances. 

(b) It is believed, and supported by studies of mantle samples, that refractory lithophile elements inside Earth are 
in the same proportion as in chondritic meteorites. In this way, primitive abundances of Th and U can be 
derived by rescaling the chondritic values. 

(c) Potassium, being a moderately volatile elements, could have escaped in the planetesimal formation phase. Its 
absolute abundance is best derived from the practically constant mass ratio with respect to uranium observed 
in crustal and mantle derived rocks. 

There are several calculations of element abundances in the BSE model, all consistent with each other to the 
level of 10%. By taking the average of results present in the literature, in ( Mantovani et al\ . 12004 ) the following 



values were adopted^^: for the uranium abundance aBSE(U) = 2 x 10 for the ratio of elemental abundances 
Th/ U = aBSE(Th)/aBSE(U) = 3.9, and K/U = aBSE(K)/aBSE(U) = 1-14 x 10"*. For a comparison, a recent re- 
view (jPalme and O'Neilll . [20031 ) - subsequent to ([Mantovani et aZ.l . l2004l ) - gives aBSE(U) = 2.18 x 10"^ Th/U = 3.83, 
and K/U = 1.2 x 10''. 



C. The crust 

Earth is the only planet, in our solar system, that has both liquid water and a topographically bimodal crust, 
consisting o f low- lying higher-density basa ltic oceanic crust (OC) and high-standing lower-density andesitic continental 
crust (CC) (" Rudnick and Fountainl . ll995Y 

Although the continental crust is insignificant in terms of mass (about half of a percent of the total Earth), it forms 
an important reservoir for many of the trace elements on our planet, including the heat producing elements. It also 
provides us with a rich geologic history: the oldest dated crustal rocks formed within 500 Ma (million years) of Earth 
accretion, whereas the oceanic crust records only the last 200 Ma of Earth history. 



We shall always refer to element abundances in mass and we remind the reader that the natural isotopic composition is ^'^^U/U = 0.993, 
232Th/Th = 1, and "Ok/K = 1.2 x 10""' 
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The crust extends vertically from the Earth's surface to the Mohorovicic (Moho) discontinuity, a jump in compres- 
sional wave speeds from « 7 to « 8 km/s which occurs, on the average, at a depth of « 40 km for the continental 
crust and at a depth of about 8 km for the oceanic crust. 

The Conrad discontinuity separ ates the continental crust in to two parts. Actually, based on additional seismic 
information several authors, e.g., (jRudnick and Fountainl . Il995l ). identify three components in the crust, the upper- 
, middle-, and lower-crustal layers (which we shall refer to as UC, MC and LC, respectively). The upper crust 
is readily accessible to sampling and robust estimates of its composition are available for most elements, whereas 
the deeper reaches of the crust are more difficult to study, so that the estimated elemental abundances are more 
uncertain. The observations show that the crust becomes more mafic"'^* with depth and the concentration of heat 
producing elements drops rapidly from the surface downwards. Not only the crust is vertically stratified in terms of 
its chemical composition, but it is also heterogeneous from place to place. This makes it difficult to determine the 
average composition of such a heterogeneous mass. 



1. Abundances of heat generating elements 

For each component of the crust, one has to adopt a value for the abundances^^ a(U), a(Th), and a(K) and to 
associate it with an uncertainty. In the literature of the last twenty years one can find many estimates of abundances 
for the various components of the crust (upper, mid dle, lower crust = UC, MC, LC and oceanic cr ust = OC), generally 
without an error value, two classical reviews being l|Tavlor an d McLennanl . [T99llWedepohil . ll995[ ). Aver age elemental 
abundances in the continental crust, and the ir vertical distribu t ion in the three main identifiable layers have been 
presented in a recent comprehensive review (jRudnick and Gaol . l2003l ). together with a wealth of data and with a 



critical survey of earlier literature on the subject. A most useful and continuously updated source is provided by the 
GERM Reservoir database. Table IIXI presents the uranium abundances used for geo-neutrino calculation in a few 
studies. 



TABLE IX Uranium mass abundances in the Earth's reservoirs and in the silicate Earth (BSE) used in recent geo-neutrino 
studies. Units are ^g/g (ppm) for the crust, ng/g (ppb) otherwise. 

Mantovani et al. Fogli et al. Enomoto 

fMantovani et al. 2004) (Fodi et al. 2005) (Enomoto. 2005) 



Reservoir 


Units 


Adopted value 






Adopted value 


0" 


Adopted value 


UC 


ppm 


2.5 


0.3 


0.13 


2.7 


0.6 


2.8 


MC 


ppm 


1.6 






1.3 


0.4 


1.6 


LC 


ppm 


0.63 


0.45 


0.23 


0.2 


0.08 


0.2 


OC 


ppm 


0.1 






0.1 


0.03 


0.1 


UM 


ppb 


6.5 


1.5 


1.5 


3.95 


1.2 


12 


BSE 


ppb 


20 


2.5 


1.0 


17.3 


4.7 


20 



Earlier paper on geo-neutrinos adopted abundances from some review papers, without tackling the problem of 
the a ssociated uncertainties. A similar approach is taken in the recent paper by Enomoto et al. ( Enomoto et al\ . 
I2OO5I ) where the values from ( Rudnick and Fountainl . [l995h are used directly, without any estimate of the associated 
uncertainties. 

Our group adopted as reference values for the abundances the average (a) of values which were available in the 
GERM database^" in 2003, considering the spread of the reported abundances (cmax — amin)/2 as indication of the 
corresponding uncertainty. In Table HXl we also presents the standard deviation of the average: 



N 



(36) 



Mafic is used for silicate minerals, magmas, and rocks which are relatively high in the heavier elements. The term is derived from 
"magnesium,, and "ferrum" (Latin for iron), but mafic magmas are also relatively enriched in calcium and sodium. 
Throughout the paper the term abundance refers to abundance in mass. 

Geochemical Earth Reference Model (GERM) available online at jhttp : //earthref . org/GERM/l 
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Fogli et al. ( Fogli "oZI . l2005l ) basically adopt the results of a recent and comprehensive review (jR udnick and Gaol . 



I2OO3I ) for the UC, MC , LC abundances of (U, Th, K) and the uncertainties quoted in that paper. Since no error 
estimates are given in (jRudnick and for the LC, Fogli et al. assume a fractional la errors of 40%. A more 



extensive comment on uncertainties can be found in Appendix [Cl 

From the table one sees that the values adopted by different groups are generally in agreement within the quoted 
uncertainties, however with the following remarks: 

• A major difference lies in the abundance for the (poorly constrained) lower portion of the crust. In the literature, 
values as low as 0.2 ppm and as high as 1.1 ppm have been reported, as a consequence of different assumptions 
for the fraction of (uranium poor) metaigneous rocks and (uranium rich) metapelitic rocks. 



Concerning the upper crust, the values quoted by different authors using different methods (surface exposure 
data, sedimentary data and loess correlations wit h La) are consistent within about 10%. From a study of loess^^ 
correlations with La, in ( Rudnick and Gaol 120031 ) a concordant average value has being obtained, however with 



a 1(7 uncertainty of 21% mainly due to the variability of the U/La correlation. 

From the table, it emerges that the contributions to geo-neutrino production from different portions of the Earth's 
crust are markedly different, the continental crust being an order of magnitude richer in heat generating material 
than the oceanic part. Relative uncertainties, as natural, increase with depth and their assessment is at the moment 
somehow tentative. 



2. The distribution of heat generating elements 

The earlier geo-neutrino studies considered the distribution of heat generating elements as spherically symmetrical 
over the Earth's crust. Actually one has to distinguish between continental and oceanic crust since they have quite 
different contents of heat generati ng elements. In additio n the thickness of the crust significantly differs from place to 
place. More recent studies, since ()Rothschild et aLl . ll998l ). ta ke into account the actual inhomogeneity of Earth crust. 

A global crustal model on a 2° x 2° degree grid, available at ( Laske et aLl . l200TI ). has been widely used in recent years. 



Data gathered from seismic experiments were averaged globally for similar geological and tectonic set tings (such as 



Archean, early Proteroz oic, rifts, etc.). The sedimentary thickness is based on the recent compilation by (jBassin et all 
I2OOOI: iLaske et aUl200ll ). Bathymetry and topography is that of ET0P05. 

Within each 2° x 2° degree tile, one distinguishes oceans and seawater, the continental crust, subdivided into three 
sub-layers (upper, middle, and lower), sediments and oceanic crust. For all these layers values of density and depth 
are given over the globe^^. 



D. The mantle: data, models and debate 

Sandwiched between Earth's crust and metallic core , the mantl e is a 2900 km layer of pressurized rock at high 
temperature. As reviewed by Hofmann in ( Hofmannl . Il997l l2003l) . mantle models can be divided into two broad 
classes, essentially corresponding to the presently contradictory geochemical and geophysical evidence of Earth's 
interior. 



1. Geochemical and geophysical evidences 



Geochemists have long insisted on a two-layer model, in which the mantle consists of a relatively primitive layer 
below a depth of about 670 km and an upper layer that is highly depleted of heat producing elements (panel (a) in 
Fig. The two layers are viewed as separate sources of the Mid-Ocean- Ridge Basalts (MORB), which come from 
mantle regions that have been already depleted in incompatible elements by extraction of the continental crust, and 
of Ocean Island Basalts (OIB), which form by melting of deeper, less depleted or even enriched mantle sectors. Also, 
a more primitive deep layer is needed from global constraints, otherwise the amount of radiogenic elements inside 
Earth is much too small with respect to that estimated within the Bulk Silicate Earth (BSE) paradigm. 



Loess is a deposit of silt (sediment with particles 2-64 microns in diameter) that have been laid down by wind action. 
Note that additional useful databases are available at http://mahi.ucsd.edu/Gabi/reni.htinl. 
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On the other hand, over the past several years seismic tomography has provided increasingly detailed images of 
apparently cold slab descending into the deep mantle, below the 670 km boundary. If cold slabs descend into the deep 
mantle, there must be a corresponding upward flow of deep-mantle material to shallow levels (panel (b) in Fig. I16p . 
If this circulation reaches the bottom of the mantle (whole mantle convection), it would destroy any compositional 
layering below the crust in a few hundred million years (at a typical speed of 3 cm yr^^ it takes about 10^ yr to move 
down to 2900 km). 

a b 

Ocean Continent Ocean Continent 

I 



B60s 



,>S;,_.^._ -ft... 





2,900 

Liquid outer core 

FIG. 16 Models of mantle circulation, adapted from (jHofmannl . l2003l l : (a) is the traditional two-layer model with demarcation 
at 670 km and nearly complete isolation between upper and lower layers; (b) is a fully mixed model. 

In brief, the composition and circulation inside Earth's mantle is the subject of a strong and so far unresolved 
debate between geochemists and geophysicists. Geochemical evidence has been used to support the existence of two 
compositionally distinct reservoirs in the mantle, the borders between them being usually placed at a depth near 
ho = 670 km, whereas geophysics presents evidence of mantle convection extending well beyond this depth. If this 
convection involves the whole mantle, it would have destroyed any pre-existing layering, in conflict with geochemical 
evidence. 

When building their respective referenc e models for geo- neutrino production, our group in i Mantovani et aZ.l . l2004 ). 



as well as Fogli et al. I Fogli et al\ . l2005f ) . used a two-reservoir mantle model. Observational values are adopted for 



the upper mantle whereas th e lower mantle abund ances are inferred from the BSE mass balance constraint. On the 
other hand, Enomoto et al. { Enomoto et al\ . l2005l ) prefer a wholly mixed mantle, with uniform abundance within it 



derived from the BSE constraint, see Table ITXl 



2. A class of two-reservoir models 



More generally, new views on mantle convection models overcome the widely diffused model of two-layer mantle 
convection, namely an outgassed and depleted upper layer overlying a deeper, relatively primordial and undegassed 
mantle layer. The ensemble of geochemical and geophysical evidence along with terrestrial heat flow- heat production 
balance argues against both whole mantle convection and layering at 670 km depth models, suggesting the existence 
of a tr ansition between the two r e servoirs (outgassed and depleted - degassed a nd primordia l) at 1600-2000 km 
depth (jAlbarede and van der Hilstl . ITool I van der Hilst and Karaso E [l999; Kellogg et all I1999D . In the numerical 



simulation of their mantle convection model, Kellogg et al. ([Kellogg et ali . il99S ) located this boundary at a depth 
of about 1600 km. 

In order to consider the implications of the prese nt debate on mantle cir culation and composition on the predicted 
geo-neutrino fluxes, our group has also considered ( Fiorentini et aZ.I . l2005c|) the uranium distribution in a wider class 



of models, including the extreme geochemical and geophysical models, in terms of just one free parameter, the depth 
h marking the borders between the two hypothetical reservoirs (see Fig. [T7)) : 

(i) above h one assumes uniform uranium abundance in the range from 2 to 7.1 ppb, as deduced from measures of 
the depleted upper mantle. 

(ii) Below h one assumes an uniform abundance, determined by requiring mass balance for the whole Earth. This 
means that uranium mass below the critical depth, myh, is obtained by subtracting from the total BSE estimated 
mass (jtibse) the quantity observationally determined in the crust (mc) and that contained in the mantle above 
h {m<h)- 
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rriyh = mesE - mc - m<h ■ (37) 

For each model, the calculated uranium abundance in the lower portion of the mantle is shown in Fig. 1181 
This class of models includes a fully mixed mantle (MIX), which is obtained for ft, = 25 km (i.e. just below a mean 
crust thickness obtained averaging the vales for continental and oceanic crust) so that the strongly impoverished 
mantle has a negligible thickness. The traditional geochemical model (REF) corresponds to h = h^. As h increases, 
the depleted region extends deeper inside the Earth and - due to mass balance - the innermost part of the mantle 
becomes richer and closer in composition to the primitive mantle. 

Among all possible models, the case h = 1630 km is particularly interesting. Below this depth the resulting uranium 
abundance is 20 ppb, corresponding to the BSE estimate. The innermost part of the mantle is thus primitive in its 
trace ele ment composition and th e crust enrichment is obtained at expenses of the mantle content above h. Again 
following ( Fiorentini et a/.l . l2005cfl . we shall refer to this model as PRIM. 



Concerning geo-neutrino fluxes from the mantle, all the models proposed above have the same amount of heat/anti- 
neutrino sources and only the geometrical distribution is varied. The largest flux corresponds to the model with 
sources closest to the surface, i.e. to the MIX model. On the other hand, the minimal prediction is obtained when 
the sources are concentrated at larger depth, which corresponds to the PRIM case^^. 



V. TERRESTRIAL HEAT 

Earth emits a tiny heat flux with an average value q « (60 90) mW/m^, definitely smaller than the radiation 
coming from the Sun, 1.4 kW/m^, larger, however, than the energy deposited by cosmic rays, about 10~^ W/m^. 
When integrated over the Earth's surface, the tiny flux translates into a huge heat flow, H « (30 45) TW, the 
equivalent of ten thousand nuclear power plants. In this section we briefly review the estimates of terrestrial heat 
flow and discuss the sources which can sustain it. 



A. Heat flow from the Earth: data and models 



A frequently quo ted value for the total heat release from Earth's surface is the estimate by Pollack et al. 
( Pollack et aLl . [T"993[ ). H = 44.2 ± 1.0 TW. In spite of the sm all quoted error, uncertaintie s seem to be much bigger: 



a recent revisitation of the problem by Hofmeister and Criss ( Hofmeister and Crissl . [2005a( ) yields 31 ± 1 TW with a 
central value close to that quoted in the seventies. 

The heat flux q is determined by using the conduction law (if one assumes that conduction is the main mechanism 
for heat transport): one measures the temperature gradient VT in near-surface rocks and their thermal conductivity 
k, and derives: 

q ^ -kVT . (38) 

In this way one obtains a mean heat flux of 65 mW/m^ from the continents. This commonly agreed value yields a 
continental contribution to terrestrial heat flow of 13 TW. Direct measurements from the oceanic crust give comparable 
values for the flux, so that, when weighted with the surface, the oceanic contribution would be about 20 TW. By 
summing these contributions one finds a total heat flow of about 30 TW, a value commonly found in the old literature. 

However, the estimate of the actual heat flow for the ocean is more difficult and controversial. The point is that 
the heat flow determined from Eq. psp is a lower bound for its actual value, since in porous and permeable rocks 
heat may be also carried out by convective flow of interstitial fluid (water). This convection lowers the temperature 
gradient below the value it would have if the rocks were dry or impervious to water. On these grounds, other attempts 
have been developed for estimating the oceanic heat flow. 

Instead of using the distribution of conductive heat measurements, Pollack et al. have based the estimate of the 
heat flux from the oceanic crust on theoretical thermal models, such as the Half Space Cooling (HSC) model and its 
variants. The models aim at a description of both ocean depths and heat flow versus age data. In the HSC model, 
as an example, depth and heat flow vary as the square root of age and the reciprocal of the square root of age, 
respectively. The models reproduce data at a semi-quantitative level; however, the predicted heat fluxes come out to 
be larger than those provided by the conductive heat measurements, particularly for very young ages. As a result. 



This is a part of the proximity argument, which we shall exploit extensively later, see Section I Villi 
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FIG. 18 Uranium a bundance in the lower part of the mantle as a function of the critical depth h from Earth's surface, 
from ijFiorentini et al.. . .20050 ) . 
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the mean heat flow from the oceans amounts to 101 mW/m^, with the oceans giving a contribution to the terrestrial 
heat flow of about 31 TW, so that the total flow is near 45 TW. 

Note that with respect to heat flow inferred from conductive heat measurements in the ocean, there is a difference of 
11 TW which is attributed to hydrothermal flow. As already m entioned, however , this procedure has been crit icized by 
Hofmeister and Criss, with their paper opening a debate, see (jvon Herzen et Ml, I2OO5I: iHofmeister and Crissl . [2005b) . 

In conclusion, it seems to us that for the global heat flow 30 TW is a sound lower limit based on direct observations, 
whereas 45 TW is a reasonable upper limit, as it corresponds to the highest estimate available in the literature. 



B. Energy sources 

Coming to the sources of heat flow, the situation is even more complex. A comparison between the Sun and 
Earth energy inventories may be useful for illustrating the differences between the two cases and for appreciating 
the difficulties when discussing Earth's energetics. Clearly, a heat flow H can be sustained for a time t provided 
that an energy source of at least U = H x t is available. For the Sun Uq = HqIq « 5 x 10*'^ J: clearly neither 
gravitation, Ug ~ GMq/Rq « 4 x lO'^^ J, nor chemical reactions, fJchom ^ 0.1 eV x Nq « 2 x 10^''' J (where Nq 
is the number of nucleons in the Sun) are sufficient; only nuclear energy, C/„uci ~ 1 MeV x Nq « 2 x 10**^ J, can 
sustain the solar luminosity over the solar age, as beautifully demonstrated b y Gallium experiments in the previous 
decade ( Abdurashitov et oil . l2002l : lAltmann et all I2OOOI : iHampel et all [T999h . On the other hand, for the Earth one 



finds [7(3 « 4 X 10'^^ J, C/chom ~ 6 x 10'^^ J, and C/nuci ~ 6 x 10'^° J, (assuming that some 10 ^ of Earth mass consists 
of radioactive nuclei), so that any of the previous mechanisms can sustain the present heat flow for the Earth's age: 
U(s = ffffitffi « 5 X 1030 J. 

In order to understand the energetics of the Earth one has to clarify the roles of the different energy sources, 
their locati ons and when th ey have been at work. In 1980, at the end of a review on the Earth energy sources 
Verhoogen (|Verhoogenl . ll980l ) summarized the situation with the following words: 



What emerges from this morass of fragmentary and uncertain data is that radioactivity by itself could 
plausibly account for at least 60 percent, if not 100 percent, of the Earth's heat output. If one adds the 
greater rate of radiogenic heat production in the past, . . . possible release of gravitational energy (original 
heat, separation of core, separation of inner core, tidal friction . . . meteoritic impact . . . ), the total supply 
of energy may seem embarrassingly large. . . . Most, if not all of the figures mentioned above are uncertain 
by a factor of at least two, so that disentangling contributions from the several sources is not an easy 
problem. 

Anderson ( lAndersonl . l2005t ) opens a recent review, entitled "Energetics of the Earth and the Missing Heat Source 
Mistery" , with the following words: 

Global heat flow estimates range from 30 to 44 TW. Estimates of the radiogenic contribution (from the 
decay of U, Th, and K in the mantle), based on cosmochemical considerations, vary from 19 to 31 TW. 
Thus, there is either a good balance between current input and output, . . .or there is a serious missing 
heat source problem, up to a deficit of 25 TW. 

Anderson summarizes in TablelXlthe potential contributions to the terrestrial energy budget. Similar to Verhoogen, 
he notes that the potential supply from radiogenic and non radiogenic sources, up to 66 TW, can even exceed the 
observed heat flow, so that paradoxes such as the missing heat source problem can be traced to non-realistic assumptions 
and initial and boundary conditions and The bottom line is that there appears to be no mismatch between observed 
heat flow and plausible sources of heating. He also notes, and we agree, that uncertainties on the different ingredients 
of the energy balance (total outflow, amounts of radiogenic material in the Earth, . . . ) are much larger than it was 
estimated in the past. 



C. Radiogenic heat: the BSE, unorthodox and even heretical Earth models 

We recall that the canonical Bulk Silicate Earth model predicts a present radiogenic heat production of 20 TW. The 
BSE is a fundamental geochemical paradigm: it is consistent with most observations, which however regard the crust 
and the uppermost portion of the mantle only. On the grounds of available geochemical and/or geophysical data, 
however, one cannot exclude that radioactivity in the present Earth provides a larger contribution to the terrestrial 
heat flow, sufficient to account for even the highest estimate of terrestrial heat flow. 

For a comparison, let us summarize some - less orthodox or even heretical - alternatives to the BSE. 
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TABLE X Sources of thermal energy in the Earth's interior, adapted from l| Andersonl . 120051 ') . 



Energy supply (potential contributions) 


1 vv 




Non radiogenic: 






Conducted from core 


8.6 




Mantle differentiation 


0.6 




Thermal contraction 


2.1 




Earthquake mduced gravitational energy 


2 




Ftadiated seismic energy 


0.3 




Tidal friction 


1 ^ 2 




Total (non radiogenic) 


15 


^ 16 


Radiogenic: 






Present radiogenic 


19^31 




Delayed radiogenic 






(1 to 2 Ga delay between production and arrival at surface) 


5 




Total (radiogenic) 


24 


^36 


Secular cooling (0 ^ 100 K/Ga) 


0- 


^ 14 


Total input 


39 


^66 



(a) It is conceivable that the original material from which the Earth formed is not exactly the same as inferred 
from Cl-chondrites. A mode l with initial composit i on as that of enstat ite chondrites could account for a present 
production of some 30 TW (|Hofmeister and Crisd . l2005at |Javovlll995f ). 



(b) A model where the BSE abundances of U, Th, and K are proportionally rescaled by a factor of 2.2 cannot be 
excluded by the observational data, if one assumes that the missing radiogenic material is hidden below the 
upper mantle. This model gives a present radiogenic heat production of 44 TW, which matches the highest 
estimate of the present heat flow. 



(c) Starting with ( Hall and Rama Murthvl . Il97ll : iLewisl . Il97lt) . several authors have been considering the possi- 



bility that a large amount of potassium is sequestered into the Earth's core, where it could provide the light 
element that accounts for the right core density, the energy source for driving the terrestrial dynamo, and - 
more gen erally - an addition a l con tribution to Earth energy budget. This possibility has been recently re- 
vived in (jRama Murthv et al\ . l2003l) . where from high-pressure and high-temperature data it was shown that 
potassium can enter iron sulphide melts in a strongly temperature-dependent fashion so that ''^"K could be as a 
substantial heat source in the core of the Earth. 



(d) Herndon (|Herndonl . |2003) has proposed that a large drop of uranium has been collected at the center of the 



Earth, forming a natural 3-6 TW breeder reactor, see also (|Raghavanl . l2002h . In this case nuclear fission should 



provide the energy source for terrestrial magnetic field, a contribution to missing heat, and the source of the 
anomalous •^He/'^He flow from Earth. 

In summary, an unambiguous and observationally based determination of the radiogenic heat production would 
provide an important contribution for understanding Earth's energetics. It requires to determine how much uranium, 
thorium and potassium are on the Earth, quantities which are strictly related to the anti-neutrino luminosities from 
these elements. 



VI. THE REFERENCE MODEL 

A reference model for geo-neutrino production is a necessary starting point for studying the potential and expec- 
tations of detectors at different locations. 

By definition, it should incorporate the best available geochemical and geophysical information on our planet. In 
practice, it has to be based on selected geophysical and geochemical data and models (when available), on plausible 
hypot heses (when possib le), and admittedly on arbitrary assumptions (when unavoidable). These duly cautious words 
from ( Fogli et al\ . l2005t ) explain the difficulties and to some extent the arbitrariness when building such models. 
They also mean that estimates of uncertainties on the predicted geo-neutrino fluxes are at least as important as the 
predicted values. 
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A. Comparison among different calculations 

Recently a few suc h models have been presented in the literature ( Enomoto et all l2005l : iFogli et all . 120051 : 
iMantovani et alL |2004[ ) . Predictions by different authors for a few locations are compared in Table IXII 



TABLE XI Predicted geo-neutrino rate from U + Th at various locations. Rates are in TNU. All calculations are normalized 
to a survival probab i lity (P ee) ~ 0.57. For Mantovani et al. the uncertainties are estimated as (A'^i^igh — A'^iow)/6, see Table XII 
of l|Mantovani et all |2004| ) . 



Location 


Mantovani et al. 
CMantovani et al. 2004) 


Fogli et al. 
fFodi et al. 2005) 


Enomoto " 


Hawaii 


12.5 ±3.6 


13.4 ±2.2 


13.4 


Kamioka 


34.8 ±5.9 


31.6 ±2.5 


36.5 


Gran Sasso 


40.5 ± 6.5 


40.5 ±2.9 


43.1 


Sudbury 


49.6 ±7.3 


47.9 ±3.2 


50.4 


Phyasalami 


52.4 ±7.6 


49.9 ± 3.4 


52.4 


Baksan 


51.9 ±7.6 


50.7 ±3.4 


55.0 



"Private communication. 



All these models rely on the geophysical 2° x 2° crustal map of ^Bassin et al. '. jgOOOl : iLaske et all l200lh and on the 
density profile of the mantle as given by PREM (Dziewonski and An derson. 198 1 ]^ 

Concerning the adopted abundances in the crust layers, Mantovan i et al. (Mantovani et a/"] . l2004[ ) use average 
values from results available in the literature in 200 2 Fog h et al (|Fogn et all |2005|L refer to the values of the 
recent review by Rudnick and Gao i Rudnick and Gaol . l2003l i. wher eas Enomoto et al. jEnomoto et I2005D adopt 
the values reported in the 1995 review by Rudnick and Fontaine ( Rudnick and Fountain Tl995|). Mantovani et al . 
and Fogli et al. assume a ch emic ally layered mantle, with abundances in the upper ru antle from (jjochum et al\ . Il983l : 
IZartman and Hainesl . ll988[) and ( Salter and Strackg . 12004 [Workman and Haral2005l ). respectively; whereas Enomoto 
et al. consider a chemically homogeneous mantle. The adopted uranium abundances in the various reservoirs are 
compared in Table [IXl All papers use the BSE mass constraint in order to determine the adopted abundances in the 
lower portion of the mantle. 

Concerning uncertai nties on the abundances in the crust and in the upper mantle, Mantovani et al. 

( Mantovani et "oZI. I2004D estimate them fron i the spread of published values, whereas Fogli et al. ijFogli et al\ . l2005l ) 
use the errors quoted by Rudnick and Gao ( Rudnick and Gaol , l2003f ). where available, see Table US The uncertain- 
ties of the abundances in the lower mantle are obtained by Fogli et al. by propagating in the mass constrain also 
the uncertainties estimated for the BSE. Mantovani et al. release the BSE constraint and take a more conservative 
attitude, including as extreme values the possibility that the lower mantle has the same small abundances observed 
in the upper mantle and, on the other side, the possibilities that the amount of heat generating elements can sustain 
a fully radiogenic 44 TW heat flow, with most of the material being hidden in the unexplored lower mantle. The 
BSE mass constraint used by Fogli et al fixes the total amount of heat generating elements to the level of ±14% 
(Icr), whereas the range considered in Mantovanti et al. is much wider. This is at the origin of the differences in the 
uncertainties quoted in Table IXII Indeed, if the calculation of Mantovani et al. is restricted to the BSE range, the 
uncertainties become comparable to those of Fogli et al. . 

One has to remark^^ that Fogli et al. also present a systematic approach to the ubiquitous issue of covariances 
in geo-neutrino analyses. In fact, for the only (at the moment) relevant case, the correlation between U and Th, 
they use the value p = 0.94, estimated by the errors on the U/Th ratio, very close to p = 1, which was used in 
(jMantovani et al\ . \2m4\ . 



B. The contribution of the various reservoirs 



The predicted signal all over the world is shown in Fig. [111 taken from ( Fiorentini et al\ . l2005c[ ) . It presents the 
Earth as it shines in geo- neutrinos. The more intense signals arise from regions with a thick continental crust, whereas 
over the oceans the signal essentially originates from the mantle. 



A more extensive comment on the treatment of uncertainties is presented in Appendix [C] 
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FIG. 19 Predic ted geo-neutrino events from uranium and thorium decay chains, normalized to 10^^ proton yr and 100% 
efficiency, from ijFiorentini et all l2005d ^ . 




10 100 1000 10000 

R [Km] 

FIG. 20 Contributed signal as a function of distance. The percentage contribution to the event yield at Kamioka originating 
from sources within a distance R is shown for the indicated values of Am^ in units of 10~^ e V'^ at fixed sin^ 29 = .863. The 
percentage contributed neutrino flux without oscillation is also shown for comparison. From l|Mantovani et a/l 120041 '). 
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The separate contributions of the different reser voirs to the produced flux of geo-neutrinos from the uranium 
decay chain are analyzed in Table IXIIi from (,Mantovani et al] , [2004j ) . At Himalaya, a site which maximizes the 



TABLE XII Uranium: masses, radiogenic heat, and predicted fluxes. Units are 10^'' kg, TW and lO" cm"^ s"\ respectively. 
The reference values, lower and upper limits are labeled a s ref, low , and high, respectively. Crust summarizes CC and OC; UM 
(LM) denotes upper (lower) mantle. Data from (Mantova ni et aZ"l [2004'l. 









Himalaya 


Gran Sasso 


Kamioka 


Hawaii 








33° N 85° E 


42° N 14° E 


36° N 137° E 


20° N 156° W 




M(U) 


H{U) 






"1>U 




Crust low 


0.206" 


1.960 


3.337 


1.913 


1.594 


0.218 


Crust ref 


0.353 


3.354 


5.710 


3.273 


2.727 


0.373 


Crust high 


0.413 


3.920 


6.674 


3.826 


3.187 


0.436 


UM low 


0.048 


0.455 


0.146 


0.146 


0.146 


0.146 


UM ref 


0.062 


0.591 


0.189 


0.189 


0.189 


0.189 


UM high 


0.077 


0.727 


0.233 


0.233 


0.233 


0.233 


LM low 


0.147 


1.399 


0.288 


0.288 


0.288 


0.288 


LM ref 


0.389 


3.695 


0.760 


0.760 


0.760 


0.760 


LM high 


1.177 


11.182 


2.299 


2.299 


2.299 


2.299 


Total low 


0.401 


3.814 


3.770 


2.346 


2.027 


0.651 


Total ref 


0.804 


7.639 


6.659 


4.222 


3.676 


1.322 


Total high 


1.666 


15.828 


9.206 


6.358 


5.720 


2.968 



"Th is value corresponds t o an uranium abundanc e in the continental crust equal to 0.91 ppm as estimated in Ref. l lTaylor and McLennanl . 
Starting from Ref. llFiorentini et aLl . [2003) we dismiss this estim ate, since it is inconsist ent with data from all other authors, and 
we use 0.3 as the lower limit for the uranium mass in the crust, see Ref. l lFiorentini et gj] , I2004D for a discussion. The values of the fluxes 
from the crust corresponding to 0.3 are 4.92, 2.84, 2.35, and 0.33 for the four locations, respectively. 



crust contribution, the prediction is $([/) = 6.7 x 10^ cm 



whereas at Hawaii, a site which minimizes the crust 



contribution, the prediction is ^(U) — 1.3 x 10^ cm^^s^^, originated mainly from the mantle. For the Kamioka mine, 
where the KamLAND detector is in operation, the predicted uranium flux is $(J7) = 3.7 x 10^ cm~^s~^. Within 
the reference model, about 3/4 of the flux is generated f rom material in the crust and the rest mainly from the 
lower mantle. At Gran Sasso laboratory, where Borexino ( Alimonti et aLl . Il998[ ) is in preparation, the prediction is 



$([/) = 4.2 X 10^ cm~^s~^, this larger flux arising from a bigger contribution of the surrounding continental crust. A 
similar calculation for Sudbury, the place which hosts the SNO detector, gives $(J7) — 4.4 x 10^ cm^-^s~^. The crust 
contribution exceeds 80%. 

T he contribution to the si gnal as a function of the distance from the detector is shown in Fig. [201 from Mantovani et 
al. ( Mantovani et al] . \2004) for the specific case of Kamioka. One notes that some 15% of the total signal originates 
from a region within 30 km from the detector, whereas half of the signal is generated within some 600 km. We remind 
the typical linear dimension of each tile in the 2° x 2° crustal map is of order 200 km, so that any information on a 
smaller scale is lost. A better geological and geochemical description of the region surrounding the detector is needed 
for a more precise estimate of the geo-neutrino signal. 



C. The effect of uncertainties of tfie oscillation parameters 

In this review we have fixed Am^ — 8.0 x 10~^ eV^ and tan^ 9 = 0.45, which gives a n asymptotical survival 
probability {Pee) — 0.57, following the best fit of Strumia and Vissani ( Strumia and Vissanil . [20051 . The same paper 



gives a 99% CL range 7.2 x lO^-'^eV^ < Am^ < 8.9 x 10~^eV^ and 0.33 < tan^ 9 < 0.61, with the corresponding range 
for the average survival probability 0.53 < (Pee) < 0.63. 

The effect of these uncertainties on the predicted signal is presented in Table IXHIl 

The predicted signal is practically unaffected by the uncertainty on Am^: when this is varied within its 99% CL 
interval the signal changes by less than one tenth of TNU. This holds for any value of the total uranium and thorium 
mass, since the precise value of Am^ only matters in the region near the detector. In addition, we observe that the 
predictions computed for the best value (Am^ ~ 8.0 x 10~^ eV^) and for the limit Am? = oo differ by +0.3 TNU. 

The uncertainty on the mixing angle is most important: at the 99% CL the relative error on the signal AS/S « 9%, 
somehow smaller (but not negligible) in comparison with the geological uncertainties. 



33 



TABLE XIII Effect of the oscillation parameters on the predicted U + Th signal at Kamioka. The relative/absolute variation 
is computed with respect to the prediction for the best fit values (Am^ = 8.0 x 10~^ eV'^ and tan'^ 6 = 0.45). 



Parameter 




Signal variation 


Am^ = 7.2 X 10"^ 


eV^ 


+0.11 TNU 


Am? = 8.9 X 10"'' 


eV^ 


-0.09 TNU 


tan^ e = 0.61 




-7.5% 


tan^ e = 0.33 




+9.3% 



VII. REFINEMENTS OF THE REFERENCE MODEL: THE REGIONAL CONTRIBUTION 

The geo-neutrino signal depends on the total amount of heat generating elements in the Earth and on the geochem- 
ical and geophysical properties of the region around the detector. For KamLAND, we estimated that about one half 
of the signal generated in the crust comes from a region within 200 km from the detector (half of the total signal is 
originated from within 600 km). This region, although containing a globally negligible amount of heat generating ele- 
ments, produces a large contribution to the signal as a consequence of its proximity to the detector. This contribution 
has to be determined on the grounds of a detailed geochemical and geophysical study of the region. 

The study of the region around Kamioka based on a detailed analysis of the six tiles depicted in Fig. [^11 including 
the possible effects of the subducting plates below the Ja pan Arc and a discussion of the contrib ution from of the 
Ja pan Sea, has been presen ted in (jFiorentini et aLl . l2005cl and in (|Enomotol . 120051: lEnomoto et'ai\ . |2005( ) . The result 



of (jFiorentini et all l2005dl ) for this regional contribution to the signal from uranium geo-neutrinos is: 

^rcg(U) (15.41 + 3.07) TNU (3cr) . (39) 

This result is obtained by including several effects, discussed in the next subsections. These refinements increase the 
signal by about 1 TNU. The global error is obtained by adding in quadrature the individual independent uncertainties. 

The results of that paper, which only considered geo-neutrinos from uranium decay chains, have been extended to 
include geo-neutrinos from thorium. 



A. The six tiles near Kamland 



The depth d i stribu tion of the Conrad and Moho discontinuities beneath the whole of the Japan Islands are derived 
m (|Zhao et a/.l . Il992l) . with an estimated standard error of ±1 km over most of Japan territory. This allows distin- 



guishing two layers in the crust: an upper crust extending down to the Conrad and a lower part down to the Moho 
discontinuity. 

In Fiorentini et al. ( Fiorentini et ^,[20053) a map of uranium abundance in the upper crust has been built, under 



the i mportant assumptio n that the composition of the whole upper crust is the same as that inferred in Togashi et 
al. ( Togashi et all l2000l ) from the study of the exposed portion, see Fig. [2TJ There is no specific study of the lower 



part of the Japan crust, however, it is well known that there are similarities between the compositi on of the Japanes e 
crust and that of the Sino-Korean block. In an extensive compositional study of East China crust ( Gao et aZl . ll998( ). 
the uranium abundance in the lower part is estimated between 0.63 and 1.08 ppm. On these grounds Fiorentini et al. 
(|_Fiorentini et at, 2005d) adopt for the abundance in the lower crust of Japan ai,c{U) — (0.85 ± 0.23) ppm. 

We remind that for building the refer ence model, the entire Earth's crust was subdivided into 2° x 2° tiles. In 
Fiorentini et al. (Fiorentini et ^,[20053) a more detailed grid was used where each tile is subdivided into 64 (1/4)° x 
(1/4)° cells, thus with a linear scale of about 20 km. This size is chosen since the sampling density for the study 
of the upper crust in the region near Kamioka is about one specimen per 400 km^ and also, concerning the vertical 
distribution of abundances in the crust, it is presently impossible to have information on the chemical composition 
on a scale smaller than the Conrad depth, generally lying at about 20 km. 

The contributions from the six tiles to the uranium signal is <S'6(U) ~ 12.74 TNU. The calculation for the same 
region within the reference model gave S'6(U) = 14.10 TNU. A similar reduction was also found in (Enomotg, .2005i ). 
The difference is understood in terms of the depletion of U concentration in the Japanese islands with respect to the 
average continental crust, already noted in Togashi et al. ( Togashi et a/.l . [2000[) . 



With respect to the prediction of the reference model, the six tiles contribute 45% of the total signal. In more 
detail, the tile hosting Kamioka generates 30% of the total produced signal. Note that the uranium mass contained 
in the six tiles is about me = 3.3 x 10"'^^ kg, really negligible (less then 0.05%) with respect to that estimated for the 
whole Earth. 
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Fiorentini et al. ( Fiorentini et al\ . l2005dl ) consider several sources of the uncertainties affecting this estimate: 
measurement errors of the chemical analysis, discretization of the upper crust, chemical composition of the lower 
crust and crustal depth. Their effects are summarized in Table IXIVI 



TABLE XIV Errors from the regional geophysical and geochemical uncertainties. 



Source 


AS (TNU) 


Remarks 


Composition of upper-crust samples 


0.96 


3(T error 


Upper-crust discretization 


1.68 




Lower-crust composition 


0.82 


Full range 


Crustal depths 


0.72 


3(T error 


Subducting slab 


2.10 


Full range 


Japan Sea 


0.31 


Full Range 


Total 


3.07 





B. Effect of the subducting slab beneath Japan 

The Japan arc, at the crossing among the Eurasian, Philippine and Pacific plates, is the theater of important 
subduction processes. The Philippine plate is moving towards the Eurasia plate at about 40 mm/yr and is subducting 
beneath the southern part of Japan. The Pacific Plate is moving in roughly the same direction at about 80 mm/yr 
and is subducting beneath the northern half of Japan. 



In o rder to estimate the effect of the subducting slab on geo-neutrino production, Fiorentini et al. ([Fiorentini et ali 
l2005dD considered two extreme cases: (a) one assumes that the slab keeps its trace elements while subducting; (b) 
at the other extreme, it is possible that, as the slab advances, all uranium from the subducting crust is dissolved in 
fluids during dehydration reactions and accumulates in the lower part of the continental crust of Japan, thus strongly 
enriching it. 

As there is no argument for deciding which of the extreme cases (a) or (b) is clo ser to reality and in or der to 
encompass both of them, the contribution from the subducting slab was estimated in ( Fiorentini et aLl . l2005dl ) as^^: 



5siab(U) = (2.3±2.I) TNU (3a) . (40) 



Enomoto et al. l|Enomoto et al\ . [200^ only consider case (a) and get a smaller correction. 
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C. The crust below the Japan Sea 



The morphology of the Japan Sea is characterized by three major basins (Japan, Yamato, and Ulleung Basins). 
The crust of the Japan basin is general l y consi dered as oceanic, whereas the nature of other basins is controversial 
and debated. Again in ( Fiorentini et aLl . l2005ct )] two extreme models are considered: 



(a) following Ref. ( Bassin et all l200"ol : iLaske et all l200lh all the basins are formed with oceanic crust, extending 



down to 7 km below 1 km of sediments. 

(b) Deeper crustal depths (up to 19 km for the Oki bank) and thicker sediments layers (up to 4 km for the Ulleung 
basin) are reported in the literature, see Table IXVl By taking these values and assigning the abundances typical 
of continental crust, one maximizes geo-neutrino production. 



TABLE XV The vertical extensions (km) of crustal layers in the Yamano basin (YB), Oki bank (OK), and Ulleung basin (UB) 
used for model (b). 





YB 


OK 


UB 


Sediments 


1.2 


0.3 


4 


Upper 


2.8 


8.7 


2 


Lower 


8.5 


10.5 


8 



In order to encompass these two extreme cases, Fiorentini et al. (|Fiorentini et ^,[20053) fix the contribution to 
the signal from the Japan Sea as: 

S'js(U) = (0.37 ±0.31) TNU (3cr) . (41) 



D. Thorium contribution and the total geo-neutrino regional signal 

By adding the above contributions, and summing in quadrature independent uncertainties one obtains equation 
P9p . The same analysis repeated for Th gives a regional contribution: 

S',cg(Th) = (3.66 ±0.68) TNU (3cr) . (42) 

Assuming a complete correlation for the errors, one has: 

S'reg(U + Th) = (19.1 ±3.8) TNU (3cr) . (43) 



VIII. BEYOND THE REFERENCE MODEL 
A. Overview 

As discussed in the preceding sections, masses of heat generating elements in the Earth are estimated on the grounds 
of cosmochemical arguments, based on the compositional similarity between Earth and carbonaceous chondrites. 
Measurements of samples from the Earth's crust imply that the crust contains about one half of this global estimate, 
whereas the mantle - which should contain the rest - is practically unexplored in this respect. A direct determination 
of the mass of heat generating elements in the globe is clearly an important test of the origins of the Earth and will 
fix the radiogenic contribution to the terrestrial heat flow, which is a presently a debated issue, see Sect. |Vl 

The geo-neutrino signal depends on the total mass of heat generating elements in the Earth and on the geochemical 
and geophysical properties of the region around the detector. The region close to the detector, although containing a 
globally negligible amount of uranium, produces a large contribution to the signal as a consequence of its proximity to 
the detector. This contribution has to be determined on the grounds of a detailed geochemical and geophysical study 
of the region, if one wants to extract from the total signal the remaining part which carries the relevant information on 
the mass of heat generating elements. Such a study of the region around Kamioka has been presented in the previous 
section. 

The contribution to the geo-neutrino signal from the rest of the world depends on the total amount of heat 
generating elements as well as on their distribution inside the Earth, since the closer is the source to the detector the 
larger is its contribution to the signal. For each value of the total mass, we shall construct distributions of abundances 
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which provide the maximal and minimal signals, under the condition that they are consistent with geochemical and 
geophysical information on the globe. 

This will bring us beyond the reference model. Essentially, we shall build models of the Earth which respect the 
observational data available, concerning the abundances in the crust and the density profile of the Earth, and we shall 
release the BSE constraint on the global amounts of heat generating elements. 

In practice, we shall keep the masses of heat generating elements in the unexplored lower mantle as free parameters. 
One can still vary the abundances along the mantle depth (which are generally believed to increase with increasing 
depth) and obtain different geo -neutrino signal for the same total amount of heat generating elements. The important 
point (IPiore ntini et al\ . l2005dl ) is the following: the assumption that the abundances are spherically symmetrical 
and non-decreasing with depth will be enough to provide rather tight constraints on the mantle contribution to the 
geo-neutrino signal. 

This has to be further combined with the contribution from the crust, which also can be maximized/minimized by 
varying the abundances in the range allowed by observational data. 

The combination of this information with the regional study allows to find the connection between the geo-ncutrino 
signal and the masses of heat generating elements in the Earth. 

This relationship will be developed in the following subsections, where we elucidate the proximity argument and 
then we combine the regional contribution for Kamioka found in the previous section with that of the rest of the 
world. 

In principle, this argument can be developed separately for uranium and thorium, i.e. connecting the geo-neutrino 
signals S'(U) and S'(Th) with the respective masses to(U) and m(Th) all over Earth. 

The main result is shown in Fig. 1221 which presents for Kamioka the connection among the uranium mass, heat 
generation, and geo-neutrino signal. 

Hr(U) [TW] 




0.5 1 1.5 2 

m[10^^kg] 

FIG. 22 The predicted signal from uranium geo-neutrinos at KamLAND, adapted from ijFiorentini et al ]. l2005dh . 

By assuming the BSE mass constraint, our prediction is 

S{V) = (29.5 ± 1.6) TNU (Ict, BSE) . (44) 
If the chondritic Th/U ratio is assumed, one finds 

5(U Th) = (32 ± 2) TNU (Icr, BSE) . (45) 
However, if the BSE constraint is released we get 



^(U + Th) = (321°) TNU (Icr, noBSE) 



(46) 
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B. The proximity argument 

The main question is to build models which, for a given total uranium mass in the Earth, m(U), provide the minimal 
and maximal signals, with the additional constraint that these models be consistent with available geochemical and 
geophysical observational data. This result can be accomplished by means of what we call the "proximity argument" : 
the minimal (maximal) contributed flux is obtained by placing heat radiogenic elements as far (close) as possible to 
the detector ( Fiorentini et all . \2004) . 



This argument can be used in several steps. 

(i) For a place on or near the continental crust, since the continental crust lies on the average closer than the mantle 
to the detector the maximal (minimal) signal is obtained by putting as much (little) radiogenic material in the 
crust, as consistent with the observed values^^. This determines the mass of uranium mass in the crust, mc(U), 
which is constrained by observational data to lie in the interval (0.3 0.4) x 10^^ kg. 

(ii) This leaves us with the problem of distributing the remaining mass, m(U) — mc(U), inside the mantle. Under 
the assumptions that the abundances in this reservoir are non-decreasing functions of the depth, the extreme 
predictions for the signal are obtained by: (a) placing heat generating elements in a thin layer at the bottom, 
or (b) distributing it with uniform abundance over the mantle. 

(iii) One can then combine the extreme cases so as to obtain lower and upper limits to the global contribution to 
the signal. 



C. The case of KamLAND 

In the case of Kamioka, after excluding the region where a separate geo-chemical and geophysical investigation has 
been performed (see Sec. IVIlT) . the proximity argument provides for the signal from the rest of the world: 

5'rw(U) = (2.25 + 14.76 x to(U)) ± (-0.55 + 2.61 x m(U)) , (47) 

where the signal is in TNU, the mass is in units of 10^^ kg and the interval within the ± sign corresponds to the full 
range of models which have been considered. 

By combining this results with the regional contribution, calculated in the previous section, S'rcg(U) = (15.41 ± 
3.07) TNU, we get the uranium geo-neutrino signal as a function of uranium mass in the Earth: 

5(U) =5o(U)±A(U) , (48) 

where: 

5'o(U) 17.66+ 14.76 X m(U) (49) 

A2(U) = (3.07)2 + (2.61 X m(U) - 0.55)^ . (50) 

This error is obtained by combining in quadrature all geochemical and geophysical uncertainties discussed in the 
preceding paragraphs. All of them have been estimated so as to cover ±3cr intervals of experimental measurements 
and total ranges of theoretical predictions. 

However, this error does not account for present uncertainties on neutrino oscillation parameters. We remind (see 
Section fVLCp that the uncertainty on the mixing angle implies a 99% CL relative error on the signal AS/S « 9%, 
which is somehow smaller (but not negligible) in comparison with the geological uncertainties. 

For the sake of discussing the potential of geo-neutrinos, we shall ignore in the following the error originating from 
uncertainties on the mixing parameter, which however should be measured more accurately. 

The expected signal from uranium geo-neutrinos at KamLAND is presented as a function of the total uranium 
mass m(U) in Fig. 1221 The predicted signal as a function of m(U) is between the two lines denoted as 5'high and S'low, 
which correspond, respectively, to S'o(U) ± A(U). 

We remark that the extremes of the band correspond to the whole range of uncertainty, which is estimated according 
to the following criteria: 



The opposite holds for places very far from the continents as the Hawaii Islands, see Section IVIII.DI 
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(i) for statistical errors we consider a ±3cr interval; 

(ii) for systematic uncertainties of geochemical and geophysical parameters we determine an interval such as to 
cover all modern estimates which we found in the literature; 

(iii) independent errors are combined in quadrature. 

We remark that the "proximity argument" , combining global mass balance with geometry, is very powerful in con- 
straining the range of fluxes: in the allowed band are enclosed all the models consistent with geological data. 

Since the minimal amount of uranium in the Earth is 0.3 x 10^^ kg (corresponding to the minimal estimate for the 
crust and the assumption of negligible amount in the mantle), we expect S'(U) to be at least 19 TNU. On the other 
hand, the maximal amount of uranium tolerated by Earth's energetics^^, 1.8 x 10^^ kg, implies S(U) not exceeding 
49 TNU. 

For the central value of the BSE model, m{\J) = 0.8 x lO^^^ kg, it was found in (jFiorentini et all . l2005dl ) 5(11) = 
29.5 ± 3.4 TNU, i.e. with an accuracy of 12% at 3a. 

We remark that estimates by different authors for the uranium mass within the BSE are all between (0.7 0.9) x 
10^^ kg. This implies that the uranium signal has to be in the interval (24.7 34.5) TNU. The measurement of 
geo-neutrinos can thus provide a direct test of an important geochemical paradigm. 

We do not expect that the next generation of experiments can collect enough statistics so as to clearly separate 
the two components (U and Th) in the signal (see Section |lT| and will have to rely on the chondritic estimate for the 
ratio for the global abundances of these two elements. On these grounds, we shall also assume this value of the ratio 
and rescale the results calculated for uranium to get the total signal as a function of the radiogenic heat from U + Th. 
In this way one gets Fig. [221 



Signal H(U+Th) from Uranium and Thorium geo-neutrinos at Kamioka 
70 I 




10 r 



r 1 

10 20 30 40 

H(U+Th) [TW] 

FIG. 23 Predictions on the combined signal S(U + Th) from uranium and thorium geo-neutrinos at Kamioka as a function of 
the radiogenic heat production rate H{\J + Th). The area between the two lines (Shigh and Siow) denotes the region allowed by 
geochemical and geophysical constraints. The shaded area denotes the region allowed by the BSE constraint. Earth energetics 
implies the signal does not exceed 62 TNU, and U and Th measured in the crust imply a signal of at least 24 TNU. 



For an uranium mass m(U) = 1.8 X 10^^ kg and relative abundances as in Eq. (|4]l, the present radioge nic heat production r ate from U, 
Th, and K decays equals the maximal estimate for the present heat flow from Earth, Hmax = 44 TW l lPollack et aUFlQQSt) . 
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D. Predictions at other locations 

In order to understand the potential of detectors at other locations one should perform a study similar to the one 
just outlined for Kamioka. In the absence of a detailed geochemical and geophysical study of the region surrounding 
the detector, we can tentatively address the problem by using the prediction for the crust as derived by the 2° x 2° 
crustal map^^. 

We resort again to the proximity argument, by making it more general, so as to cover also locations which lie far 
from the continents. 



a. Contribution from the crust. The contribution, as a function of the uranium mass contained in the crust toc(U), 
can be obtained by rescaling the predictions of the reference model. This gives at the i-th location: 

Sh^a^x mc(U) , (51) 

where the response coefficients ai are presented in Table IXVll for a few locations. 



TABLE XVI The crust response coefficients for a few locations. 



Site 


a 


Hawaii 


8 


Kamioka 


57 


Gran Sasso 


69 


Sudbury 


93 



b. Contribution from the mantle. We remind that, under the assumptions that the abundances in this reservoir are 
radial and non decreasing function of the depth, the extreme predictions for the signal are obtained by: 

(i) placing uranium in a thin layer at the bottom; 

(ii) distributing it with uniform abundance over the mantle. 

For a mass tom(U) in the mantle, the two cases give, respectively, and independently of the location: 

5mjow(U) = Aow X mM(U) TNU , (52) 
5M,high(U) = /3high X mM(U) TNU , (53) 

with /3iow = 12.15 TNU and /3high — 17.37 TNU when the mass, as here and in the following, is measured in units of 
10^^ kg. 



c. Combining mantle and crust. At each location, the total uranium signal will be: 

S'(U) = S'c(U) + S'M(U) = a, xmc(U)+/3xmM(U) = (a, -/3)xmc(U) + /?xm(U) , (54) 

where the total uranium mass is m(U) — r7ic(U) + mM(U). 

If one wants the largest signal, obviously one has to put /3 = /3high and, if ai > /?high, mc(U) has to be as large 
as possible, consistently with the observations of the crust, i.e. mc(U) = min(m(U), 0.41). In the opposite case, 
c^i < /3high, one takes instead mc(U) as small as possible, i.e. mc(U) = 0.3. Note that the first case is what occurs 
on or close to the continental crust, whereas the second case corresponds to locations far from the continents. 

Similarly, for minimizing the signal one puts /3 = /3iow and then one takes mc(U) = 0.3, if ai > /Siowj or mc(U) — 
min(TO(U), 0.41) in the opposite case. 



We remind that for Kamioka the more detailed calculation changes the signal by just about 1 TNU; however the uncertainty of the 
regional contribution is about ±3 TNU (at 3cr). 
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The results for the extreme cases, shown in Fig. [24l deserve the foUowing comments. The band is most narrow for 
the Hawaii, as natural since this place, far from the continents, is most sensitive to the amount of radioactivity hidden 
in the mantle. On the other hand, it is more wide at Sudbury, where the signal is dominated from the contribution 
of the crust. Concerning Kamioka, the contour defined by this analysis are close (within four TNU or less) to those 
which were derived by adding the geological and geochemical information on the region, see Fig. [22l This gives us 
some confidence about the predictions for the other sites. 



IX. KAMLAND RESULTS AND THEIR INTERPRETATION 

A. Overview 

KamLAND (Kamioka Liquid scintillator Anti-Neutrino Detector) is so far the largest low-energy antineutrino 
detector ever built and studies a wide range of science, spanning particle physics, geophysics and astrophysics^^. 

The KamLAND collaboration is made up of researchers from Japan, China, France and United States. The detector 
is situated in the same cavern used by the original Kamiokande experiment, where 2002 Nobel laureate Masatoshi 
Koshiba performed much of his ground-breaking research in neutrino physics. 

Located on the island of Honshu in Japan, since 2002 KamLAND detects hundreds of anti-neutrinos per year 
from nuclear reactors hundreds of kilometers away, an enormous improvement over previous attempts with any 
other detector. KamLAND has observed an anti-neutrino d eficit as well as energ y spectral distort ion confirming 
neutrino oscillations and hence non vanishing neutrino masses (|Araki et all [2005b: E guchi et adl2003t) . As a natural 



continuation of the scientific program, KamLAND aims now at the direct observation of ^Be solar neutrinos by 
detecting recoil energy in neutrino-electron scattering processes. 

KamLAN D has also pe rformed a geo-neutrino investigation. In 2005, the KamLAND collaboration has pub- 
lished (jAraki et al\ . l2005at l experimental results, claiming some 28 geo-neutrino events from uranium and thorium 
decay chains in a two-year exposure. This important step shows that the technique for exploiting geo-neutrinos in 
the investigation of the Earth's interior is now available. From the KamL AND data, including ne w measurements of 
the i3C(a,n)i^0 cross section, one finds 5(U-|- Th) = (68^35) TNU, see ([Fi orentini et all l2005a[) . The central value 
is close to the prediction of a maximal and fully radiogenic model (see Sect. IVHLCp . however the BSE prediction is 
within 1(7 from it. 

In the future, with more statistics KamLAND should be capable of providing a three-sigma evidence of geo- 
neutrinos, but discrimination between BSE and fully radiogenic models definitely requires new detectors, with class 
and size similar to that of KamLAND, far away from nuclear power plants. 

In the next subsections we shall present KamLAND results on geo-neutrinos and discuss their implications on the 
terrestrial heat. 



B. The KamLAND detector 

Neutrinos are detected at KamLAND^° by the inverse beta-decay reaction, 

De+p-^ e+ +n - 1.806 MeV , (55) 

with a large amount of organic liquid scintillator. The liquid scintillator essentially consists of hydrocarbons (C„iJ2n) 
which provide the hydrogen nuclei acting as the target for antineutrinos. We remind that the energy threshold of the 
reaction, 1.806 MeV, is low enough to detect a part of the U-series and Th-series geo-neutrinos. 

The reaction produces two correlated signals. A prompt signal is given by the slowing down positron and by the two 
0.51 MeV gamma rays from positron annihilation. The delayed signal consists of a 2.2 MeV gamma particle, which 
is emitted in the thermal neutron capture on proton. The thcrmalization and capture processes take place in about 
200 microseconds, and neutron capture occurs typically 30-50 cm apart from the neutrino reaction vertex. The time 
and space correlations of the two signals are distinguishing characteristics of electron-type antineutrino events. The 
delayed coincidence of the two signals provides an effective way to select antineutrino events with excellent separation 
of background events. 



An overview of KamLAND science plan is found in l|Grattal . Il999l ). 

An extensive description of the detector can be found in the PhD thesis of S. Enomoto llEnomotol . l2005l) . 




FIG. 24 The predicted signal from uramium geo neutrinos at Hawaii (top- left), Kamioka (top-right), Gram Sasso (bottom- left), 
and Sudbury (bottom-right). The area between the black line (Shigh) and the grey line (Slow) denotes the region allowed by 
geochemical and geophysical constraints. 
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The KamLAND detector basically consists of 1000 tons of ultra-pure liquid scintillator (LS) contained in a 6.5 m 
radius spherical balloon and of 1879 surrounding photomultiplicr tubes (PMT) that cover 34% of the sphere. The 
detector is located 1000 m underground in the Kamioka mine, just beneath the Mt. Ikenoyama summit, Gifu, Japan 
(36.42° N, 137.31° E). The 2700 m water equivalent thickness of rock covering the detector sufficiently reduces cosmic 
muon flux, resulting in 0.34 Hz of muon event rate. 

Figure [25l illustrates the KamLAND detector. The LS, balloon and PMT's are contained in a 9 m radius spherical 
stainless steel vessel. 1325 17-inch PMT and 554 20-inch PMT's are mounted inside the stainless steel vessel viewing 
the center of the LS sphere. The 6.5 m radius LS-containing balloon is positioned at the center of the stainless steel 
vessel, being supported and constrained by a network of Kevlar ropes. Non-scintillating mineral oil (MO) is filled 
between the stainless steel vessel and the LS containing balloon, providing gravity/buoyancy balance to the LS sphere, 
and also acting as a buffer layer against radiations into the LS from the stainless steel vessel, PMT, and everything 
surrounding the vessel. The MO layer is further divided into two spherical shells by 8.25 m radius transparent acrylic 
wall, to isolate the balloon contacting MO from the PMT/vessel exposed MO and reduce radioactive contamination 
around the LS. 




FIG. 25 An overview of the KamLAND detector from (jEnomotd . l2005l ;i. 



The inner part of the 9 m radius stainless vessel is called the inner detector (ID), whereas the portion outside of 
the 9 m radius vessel is called the outer detector (OD). The space between the vessel and the cave wall is filled with 
3200 tons of pure water, and viewed by 225 20 inch PMT's. Cosmic muons passing through the OD are tagged by 
the OD PMT's, by detecting Cherenkov light. 



C. KamLAND results on geo-neutrinos 

The KamLAND coUaboration has reported (jAraki et aLl . l2005al ) data from an exposure of Np = (0.346±0.017) x 10^^ 
free protons over a time T = 749 days with a detection efficiency e = 68.7%; the effective exposure is thus a = 
NpXTxs^ (0.487 ± 0.025) x 10^^ proton-yr. 

In the energy region where geo-neutrinos are expected, see Fig. 1261 there are C = 152 counts, implying a statistical 
fluctuation'^^ of ±12.5. Of these counts, a number R = 80.4 ± 7.2 are attributed to reactor events, based on an 
independent analysis of higher energy data. Fake geo-neutrino events"^^, originating from -'^^C(q!, n)^^0 reactions 
following the alpha decay of contaminant ^^"Po, are estimated to be F = 42 ± 11, where the error is due to a 20% 



In this section and in the following ones the quoted errors correspond to Icr interval, 
bee Section IDCP] for a more detailed discussion of this point. 



43 



uncertainty on the ^'^C(a,n)^^0 cross section and a 14% uncertainty on the number of ^^"Po decays in the detector. 
Other minor backgrounds account for B = 4.6 ± 0.2 events. 

A straight estimate by subtraction, A^(U + Th) — C ~ R — F — B, with an uncertainty obtained by combining 
quadraticaUy errors gives: A^(U + Th) = 25 ± 18 . 

KamLAND "rate only" analysis (jAraki et aZ.I . l200"5ah . which includes detection systematic errors, partially corre- 
lated with background errors, gives A^(U + Th) = 25^28! corresponding geo-neutrino signal is thus S{\J + Th) = 
iV(U + Th)/a = 5i;i^^ TNU. 

This "rate only" study has been improved in Ref. ( Araki et all l200"5al ) by exploiting the shape of the spectrum, 
with the ratio of events A^(U)/iV(Th) being fixed at the chondritic (BSE) prediction. A likelihood analysis of the 
unbinned spectrum, see Fig. [#1 yields iV(U + Th) = 28tlt which implies S'(U + Th) 57^11 TNU. 

As a curiosity, an analysis where both N{IJ) and A^(Th) are left unconstrained yields as a best fit A^(Th)/A^(U) w 
5.7, which looks far from the chondritic value ~ 0.25; however, both values are comfortably consistent with the data 
already at the la level. In fact, the statistics is so poor that one cannot say anything from the KamLAND data 
concerning iV(Th)/iV(U), see Fig.^H 



D. Fake antineutrinos and a refinement of the analysis 

As a result of ^^^Rn contamination, ^^"Pb, which has half-life of 22 yr, is distributed throughout the detector. It 
produces (see Fig. [2]) ^^°Po which decays emitting a particles with a kinetic energy of 5.3 MeV. These act as a neutron 
source through ^^C(q!, n)^^0 reactions occurring on the "'^'^C nuclei which are present in the organic scintillator. The 
neutrons in the ^^C(a,n)^^0 reaction are produced with kinetic energy up to 7.3 MeV. Owing to scintillation-light 
quenching for high-ionization density, only about one-third of this energy is converted into visible energy as the 
neutrons thcrmalizc. The thermal neutrons are captured by protons with a mean capture time of 200 microseconds, 
producing a delayed signal identical to that from neutron inverse /3-decay. In summary, one has a fake antincutrino 
signal. 

In order to extract the true geo-neutrino signal one has to subtract these events. As already remarked, a major 
uncertainty originates from the ^^C{a,n)^^0 cross section^'^. 

The number of ^^C nuclei in the fiducial volume is determined from the measured ■^^C/^'^C ratio in the KamLAND 
scintillator. On the basis of the ^^C(a, n)^^0 reaction cross-section, the a-particle energy loss in the scintillator, 
and the number of ^^"Po decays, the total number of neutrons produced is expected to be 93 ± 22. This error is 
do minated by the unce rtainty in the total ^'^C(q!, n)^ ^0 reaction cross sec tion. The values for the cross section used 
in ( Araki et all . l2005al ) are taken from the JENDL ( Shibata et all |2002| ) compilation , which provides an R- matrix 



fit of relatively old data. A 20% overall uncertainty has been ad opted in ( Araki et oZj, l2005a[) . corresponding to the 



accuracy claimed in the original experimental papers (see, e.g., (jShibata et a l. 20021)). 

Re cently a series of high-precision measurements for this cross section has been performed ( Harissopulos et all . 
I2005D . In the relevant energy range (1 5.3) MeV, the absolute normalization has been determined with a 4% 
accuracy. The measured values are genera l ly in ve ry good agreement with those recommended in JENDL, see Fig. 1291 
however, it was found in ( Fiorentini et aLl . l2005a[) that the neutron yield per a particle is 5% smaller. It follows that 
the number of fake geo-neutrinos is lower, F = 40 ± 5.8, and geo- neutrino events obviously increase'^^ 



The "rate only" analysis gives now 27^15 geo-neutrino events ( Fiorentini et all l2005a ). corresponding to 5'(U 



Th) = 55^3^ TNU. An analysis of the binned spectrum has also been performed in ( Fiorentini et a/.l . l2005a[) with the 
result N{\J + Th) = 311^^ counts, corresponding to ^(U + Th) = 63l^^ TNU. 

These signals should be compared to S{V + Th) = 51^11 TNU and S{\J + Th) = 57tll TNU, respectively, which 
were obtained using the JENDL ^^C(a,rt)^^0 cross section. 

In summary, by using the new high-precision data on ^^C{a,n)^^0 one extracts from KamLAND data a larger 
geo-neutrino signal with a smaller error. This corroborates the evidence for geo-neutrinos in KamLAND data, which 
becomes close to 2.5(J. 



In fact, the claim of 9 geo-neutrino events in llEguchi et a;.Ll2003l) should be dismissed: more than half of these events are to be considered 
as fake sign als, produced from t he ^^C(o,n)^®0 reaction. 

Indeed Ref. llAraki et aiy\2005ei) mentions that an alternative analysis including the time structure of the scintillation light from different 
particles produced a slightly larger geo-neutrino signal, which is consistent with the result presented here. 



44 




FIG. 26 The energy spectra in KamLAND, adapted from l|Araki et all l2005al V Main panel, experimental points together 
with the total expectation (thin dotted black line). Also shown are the total expected spectrum excluding the geo- neutrino 
signal(thick solid black line), the expected signals from ^'^^U (dot-dashed line labeled U) and ^^^Th (dotted line labeled Th) 
geo-neutrinos, and the backgrounds due to reactor antineutrinos (dashed line labeled Reactor), ^'^C{a,n)^^0 reactions (dotted 
line labeled ^^C{a,n)), and random coincidences (dashed line labeled Rnd). 




FIG. 27 Confidence intervals for the number of geo-neutrinos detected, from (jAraki et al ]. l2005al '). assuming the mass ratio 
Th/U = 3.9. The shaded area represents the prediction of the reference model. 
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FIG. 28 Confidence intervals for tiie number of geo-neutrinos detected, from (|Araki et al. I l2005al ). wlien tlie mass ratio Th/U 
is left as free parameter. The shaded area represents the prediction of the reference model. 
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FIG. 29 Cross section of ^^C(a , n)^^0. The solid line corresponds to the JENDL compilation, dots are the experimental points 
from (|Harissopulos et oZl . 120051 ) . 



E. Implications of KamLAND results 

The geo-neutrino signal reported by KamLAND, ^(U + Th) = 57^3^ TNU, is well consistent with the BSE 
prediction, k, 37 TNU, as well as with the ~ 56 TNU prediction of models for maximal and fully radiogenic heat flow, 
see Fig. [23l 

In order to extract some more quantitative information from the data, we have to extend Fig. 1231 including models 
which produce even larger heat and signal. These models have been built so that an arbitrary amount of uranium and 
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thorium in the chondritic proportion is hidden in the mantle'^^. The allowed band in Fig.[30]is estimated by considering 
rather extreme models for the distributions of radioac tive elements, chosen to maximize or minimize the signal for 
a given heat production rate (jFiorentini et a/l l2005at l. We also remark that, in comparison with the experimental 
error, the width of the al lowed band is so narrow that we can limit the discussion to the median line in Fig. 1301 which 
represents, according to ( Fiorentini et a/.l . l2005al) . the best estimate for the relationship between signal and power. 
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FIG. 30 The combined signal from uranium and thorium geo-neutrinos and the radiogenic heat production rate. The shaded 
area denotes the region allowed by geochemical and geophysical constraints. The dashed median line represents our best 
estimate for the relationship between signal and power. 



By considering only this median line in Fig. [301 the reported signal S'(U+Th) = 57+3^ TNU imphes^*^ H{U + Th) = 
38133 TW (rate + spectrum) and the 99% confidence limit on the signal (145 TNU) corresponds to 133 TW. By 
including the uncertainty band of the theoretical models, one gets an upper bound of 162 TW, see point A in 
Fig. 1301 This point corresponds to a model with a total uranium mass m(U) = 8 x 10^^ kg, an uranium poor crust, 
r7ic(U) = 0.3 X 10^^ kg, the rest of the uranium being placed at the bottom of the mantle, and global chondritic 
thorium-to-uranium ratio. 

This 162 TW upper bound is much higher than the 60 TW upper bound claimed in l Araki et aZ.I . l2005al) . which 
was obtained by using a family of geological models whi ch are too narrow and a re also incompatible with well-known 
geochemical and geophysical data, see the discussion in f Fiorentini et ^,|2005a'). 

We remark that the bound H{\J + Th) < 162 TW does not add any significant information on Earth's interior, 
since anything exceeding H{IJ + Th) = 37 TW {i.e. H{\J + Th+ K) — 44 TW) is unrealistic. The upper limit simply 
reflects the large uncertainty in this pioneering experiment. 

In summary, KamLAND has shown that the technology for geo-neutrino detection is now available; however, the 
determination of radiogenic heat power from geo-neutrino measurements is still affected by a 70% uncertainty. 

An important quantity for deciding the potential of future experiments is the relationship between geo-neutrino 
signal and heat production. The basic parameter is the slope, dS/dH, which expresses how the experimental error 
translates into an uncertainty on the deduced heat production. For our models we find from Fig. dS/dH w 
1 TNU/TW. Discrimination between BSE, H{V + Th) « 16 TW, and fully radiogenic models H{V + Th) w 37 TW, 



We note that models with H{IJ + Th) > 37 TW are essentially unrealistic; this portion of the graph is included just for discussing 
KamLAND results. 

By using the result from the analysis including the new values of ^'^C{a,n)'^''0, S(U + Th) = 63^11 TNU, one finds H(V + Th) = 
44t?J, TW. 
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requires a precision AiJ w 7 TW, and thus an experiment with an accuracy AS* « 7 TNU. 



X. BACKGROUND FROM REACTOR ANTINEUTRINOS 



As first pointed out by Lagage ( Lagagel Il985( ). antincutrinos form nuclear power plants can be a significant back- 
ground for geo-neutrino detection. 

An order of magnitude estimate of the flux of antincutrinos from reactor can be immediately found from the 
knowledge of the energy produced per fission {Ens ~ 200 MeV) and the number of antincutrinos resulting form each 
fission {Nc « 6). The flux at a detector lying a distance R from a reactor with thermal power W is thus: 

For a typical value W ~ 3 GW at i? « 100 km one has 'P'^''"^ « 2.5 X 10^ cm s . In a region where there many 
reactors, as near Kamioka (21 nuclear reactors already within 200 km) the signal of man-made antincutrinos exceeds 
that from natural radioactivity in the Earth. 

In more detail, the four isotopes whose fission is the source of virtually all the reactor power are ^'^^U, ^"^^U, 
^^^Pu, and ^"^^Pu. Each isotope produces a unique neutrino spectrum through the decay of its fission fragments 
and their daughters. The instantaneous fission rates of the four isotopes are used as an input for the evaluation 
of the antineutrino spectrum. For all b ut ^^^U, car e ful measuremen ts of the (e l ectron ) spectrum from fission by 
thermal neutrons have been performed ( Hahn et all . Il989 : Schreckenbach et all Il985l ). In Fig. [31] we show the 



differential antineutrino spectrum calculated assuming 10^° fissions per second, corresponding to about 3 GW, in a 
reactor 100 km from the detector. The spectrum was calculated assuming a fuel compositio n (0.5 68,0.297,0.078,0.057) 
for (^^^■'^U, 239p^ ^ 238y^ 24ip^y . foj. 239p^ 24ip^ ^^^^ spcctra from (Hahn_eia|], I1989D . for ^ssy spectrum 
from (jSchreckenbach et all [l985l ). and for ^38^ the spectrum from IVogel et all Il98lr) corresponding to 0.5 MeV 
neutrons. For an extensive review of reactor antincutrinos see ( Bemporad et all 20021 ). 
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FIG. 31 Antineutrinos from reactors. The differential produced flux (no oscillations) in units of 2500 s~^ cm~^ MeV~^, 
inverse beta decay cross section in units of 10~*''cm'^, and the corresponding signal in TNU MeV"'^. The flux corresponds to 
10^° fissions per second (or about 3 GW) in a reactor 100 km from the detector. 

About 450 reactors are operational all over the world. If they all work at full power, this results in a total heat 
production of about 1 TW, just a factor 30 smaller than the natural heat flow from the Earth. The man-made 
antineutrino Imninosity of the Earth is thus « 10^"^ s^^, a factor 10 below the natural luminosity in gco-neutrinos 
from U and Th chains. 

Maps of the nuclear power plants in the world (see Fig. l32p . together with information on power and operational sta- 
tus are provided by several organizations, e.g., the International Nuclear Safety Center, United Nations Environment 
Programme. 
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FIG. 32 Nuclear power plants in the world, from http://www.insc.anl.gov/pwrmaps/map/world_map.php. 




FIG. 33 The ratio of reactor anti- neutrino events (in the geo-neutrino energy region) to the expected geo-neutrino events all 
over the globe. 
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The ratio r of reactor event rate in the geo-neutrino energy window to the geo-neutrino signal S'(U+ Th) predicted 
by the reference model is shown in Fig. [321 all over the globe and its inverse s = 1/r in Table IXVlII for eight locations. 
Kamioka (s w 0.14) is clearly one of the worst locations over the globe. At Gran Sasso and Sudbury the geo-neutrino 
event are comparable to reactor events, whereas a place like Hawaii looks much more favorable. 

XI. FUTURE PROSPECTS 

In summary, KamLAND has shown that the technique for exploiting geo-neutrinos in the investigation of the 
Earth's interior is now available. On the other hand, the determination of radiogenic heat power from geo-neutrino 
measurements is still affected by a 70% uncertainty. The best fit of the KamLAND result implies a radiogenic heat 
production close to the prediction of maximal and fully radiogenic model; however, the BSE prediction is within 
la from it. In order to discriminate among different models of heat production in the Earth an accuracy of at least 
±7 TW is necessary. The relationship between geo-neutrino signal and radiogenic heat, dS/dH 1 TNU/TW, implies 
that the experimental error has to be ±7 TNU, i.e. a factor of four improvement with respect to present. 

It looks to us that the following questions are relevant for the future: 

• How shall we have definite (at least Scr) evidence of geo-neutrinos? 

• How much uranium and thorium are in the Earth's crust? 

• How much in the mantle? 

• What can be said about the core? 

A preliminary point for establishing suitable detector locations is the reactor background. Figure [33] shows the 
ratio of reactor events (in the geo-neutrino energy region) to the expected geo-neutrino events all over the globe. 
KamLAND location is obviously one of the worst for the study of geo-neutrinos. 

The potential of different locations is summarized in Table IXVIIl where we present the separate contributions to 
the geo-neutrino signal rate from crust and mantle according to our reference model, S*, together with the reactor 
event rate in the geo-neutrino energy window, R. In the same Table we present two merit figures: 

(i) AS-Q = v^^Tl? is the square root of the total counts expected in a detector with an effective exposure of 
OfQ = 10'^^ proton-yr; it represents the limiting (i.e. neglecting backgrounds other than reactors, uncertainties 
of oscillation parameters . . . ) statistical error on the geo-neutrino signal which might be achieved with such a 
detector. For an exposure a the statistical error is AS* = AS'o x a^/a. 

(ii) The ratio of geo-neutrino events to reactors events in the geo-neutrino energy window s = S/R. 



TABLE XVII The geo-neutrino (U-f Th) signal rate expected from the crust Sc, from the mantle Sm, and their sum S, 
together with the reactor event rate R in the geo-neutrino energy window. ASo represents the limiting statistical error for an 
effective exposure ag — 10^^ proton-yr. All rates are in TNU. The s factor is the ratio between the geo-neutrino events and 
reactor events in the geo-neutrinos energy window. 



Location 




Sm 


S 


R 


ASo 


s 


Pyhasalmi 


42.5 


9.0 


51.5 


27.2 


8.9 


1.9 


Homestake 


42.3 


9.0 


51.3 


9.4 


7.8 


5.5 


Baksan 


41.8 


9.0 


50.8 


11.8 


7.9 


4.3 


Sudbury 


41.8 


9.0 


50.8 


53.7 


10 


0.9 


Gran Sasso 


31.7 


9.0 


40.7 


35.1 


8.7 


1.2 


Kamioka 


25.5 


9.0 


34.5 


230 


16 


0.15 


Curacao 


23.5 


9.0 


32.5 


3.3 


6.0 


9.9 


Hawaii 


3.5 


9.0 


12.5 


1.4 


3.7 


9.2 



With additional statistics KamLAND should be capable of providing 3cr evidence of geo-neutrinos, but discrimina- 
tion between BSE and fully radiogenic models definitely requires new detectors, with class and size similar to that of 
KamLAND, far away from nuclear power plants. Borexino at Gran Sasso should reach the 3a evidence, but cannot 
go much further due to its relatively small size. 
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At Sudbury, SNO+ with liquid scintillator will have excellent opportunities to determine the uranium mass in the 
crust, which accounts for about 80% of the geo-neutrino signal. This will provide an important test about models for 

the Earth's crust. 

A detector at Hawaii, very far from the continental crust, will be mainly sensitive to the mantle composition. We 
remind that the amount of radioactive materials in this reservoir is the main uncertainty of geological models of the 
Earth. Due to the absence of nearby reactors, the gco-ncutrino signal can be measured with a small error, such 
that different models for terrestrial heat generation can be discriminated. On the other hand it is necessary that 
non-reactor backgrounds be kept at very small value. 

For the very long term future, one can speculate about completely new detectors, capable of providing (moderately) 
directional information. These should allow the identification of the different geo-neutrino sources (crust, mantle and 
possibly core) in the Earth; in summary, se son rose fioriranno^^ . 
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APPENDIX A: Analytical estimates of the geo-neutrino flux 
1. The flux from a spherical shell 



Assuming spherical symmetry, the (produced) antineutrino flux at a detector on the Earth surface originated from 
a spherical shell centered at the Earth center and with radii Ri and i?2 (see Fig. [34]) is 




FIG. 34 Relevant variables for parameterizing source positions in the Earth relative to the detector D and spherical shell whose 
flux is given by Eq. ^KT\. 

From Eq. (|Aip one derives simple expressions for the contributions of the crust and of the mantle, when treating 
them as spherical shells of uniform density of heat generating elements. 
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2. Flux from the crust 



In this case i?2 — R(s and i?i = i?® — A, where A « 30 km is the thickness of the crust. Since A <C i?e, Eq. (|A1[) 
simplifies to: 



1 + log 



2Ra 



In terms of the mass of the element X contained in the crust mc{X): 

nxmc{X) 



A 



X 



mxTx^TTR'iA 



(A2) 



(A3) 



where mx and tx are the mass and lifetime of the nucleus X and nx is the number of antineutrino produced in the 
decay chain. This gives 



^c{X) 



nxmc{X) 



mxTxSTrRi 



1 + log 



2R^ 
A 



Note that this result is weakly dependent on A. 
By inserting the appropriate constants one finds 



$c(U) = 5.1 X 10'' cm^^s^^ X mc(U) , 



and 



(A4) 

(A5) 
(A6) 

(A7) 

This provides an order of magnitude estimate of the flux originated from the crust, however, as shown in Table IXlIl 
there can be substantial differences among different locations. 



$c(Th) = 1.1 X 10^ cni^^s"^ x TOc(Th) , 



where U and Th masses are measured in units of 10^^ kg. 

For the values used in the reference model (mc(U) — 0.353, mc(Th) = 1.38) one finds 



^^•(U) = 1.8 X 10^ cm^^s 



2„-l 



$c(Th) = 1.5 X 10'' cm-^s 



2„-l 



3. Flux from the mantle 

In this case R2 — i?® so that Eq. (|Aip simplifies to: 

^m{X) 
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Rm — R 
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M 



R® 2 \ i?0 — Rm , 

where Rm ~ 3500 km is the inner radius of the mantle. The specific activity is in this case 

'SnxmM{X) 



Ax = 



mxTx47r(i?| - i?3^) 



This gives 
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log 
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mxTxSnRl, + RmR® + R 
By inserting the appropriate constants one finds 

$A/(U) = 2.30 X 10^ cm-^g-i x mM(U) 

and 

$M(Th) = 0.50 X 10*^ cm^^s"^ x mM(Th) . 
For the values used in the reference model (mj\/(U) = 0.451, rriAf (Th) = 1.76) one has 

$A^(U) 1.04 X 10^ cm"2g-i $A/(Th) ^ 0.88 x 10*^ cm^^g-i 



M 



(A8) 
(A9) 

(AlO) 

(All) 
(A12) 
(A13) 



These values are in agreement with the numerical calculation used for the reference model, to the level of about 9%. 
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APPENDIX B: The contributed flux as function of the distance 

Again assuming spherical symmetry, the contribution to the flux from the portion of the crust at a distance x from 
the detector is 

d^^{ A/2 x<A 

dx I {A/2) X (A/a;) A < x < 2i?e ^ ' 

By using Eq. (|A2p the relative contribution is 

1 _ / (1 + log(2i?e/A))-i a; < A 



$c dx 1(1 + log(2i?e/A))"^ x (A/x) A < x < 2i?a 



(B2) 



This analytical estimate is shown in Fig. [35] together with the corresponding numerical result for the calculation of 
the reference model. 

A similar calculation can be made for the mantle. The result for a single- or double-shell mantle model is presented 
in Fig. [36l together with the numerical result of the reference model, which considered 25 distinct shells. 

While the uniform model overestimates the contribution from distances smaller than about 2000 km, the two-shell 
model is quite close to the reference modeP^. 



APPENDIX C: A comment on geological uncertainties 

An assessment of uncertainties is most important for understanding the significance of the theoretical predictions 
and the relevance of geo-neutrino experiments. 

In the case of geological measurements, error determination is admittedly more complex than for laboratory mea- 
surements since the quantity to be determined is often indirectly measured or extrapolated from an incomplete set of 
samples, important examples being the elemental abundances in the different Earth's reservoirs, and individual results 
are often published without quoting an error. Nevertheless it is important to have an - even if rough - estimate of 
the uncertainties and to propagate it onto the predicted signals. 

In the following we suggest an approach for estimating and combining errors of the geological quantities relevant 
for geo-neutrino calculations. 



1. Elemental abundances: selection and treatment of data. 

As for any experimental quantity, it is not possible to give completely objective criteria for the choice of data. It is 
the somewhat subjective judgment of experts in the field that selects the relevant data and uses them to extract the 
"best" educated estimate and its error. 

One has the choice of making his own compilatio n and selection of data, as in (jMantovani et "all. I2004D or of using 
some existing compilation, as in (|^ginral|,l200i). 



We believe that a robust procedure is to select all published results after excluding measurements that have been 
superseded or included in later results, or that are dependent on measurements already included, measurements that 
are clearly inconsistent with known more reliable information or that are based on questionable assumptions. Then 
all the selected measurements are averaged and the estimated standard deviation of the mean (the standard deviation 
of the results divided by \/n — 1 where n is the number of independent results) is used as error. This procedure 
implies that new independent measurements consistent with the existing ones reduce the error and new information 
can exclude inconsistent results. 

In fact this is basically the procedure adopted by the Review of Particle Physics ( Yao et all . l2006l ). with the 



addition that selected data are weighted using their errors. However, geological abundances are often published 
without standard errors: in fact it is not easy to give a reliable estimate of the error, since the final number depends 



In fact the density of uranium in the mantle in our reference model is not constant for two reason: the total density grows with 
dep th from p = 3.38 g cm~^ at R = 6346 km (depth 32 km) to p = 5.41 g cm"'' at R = 2900 km (we used 25 different lay- 
ers l|Dziewonski and Anderson) . 1 198ll )) and the abundance of uranium changes from a = 6.5 ppb in the upper mantle (above a depth of 
about 632 km) to a = 13.2 ppb in the lower mantle. 
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FIG. 35 Antineutrino flux from the crust as function of tlie distance from a detector at tfie Gran Sasso underg round laboratory, 
normalized to an unitary total flux. Points corresponds to a detailed reference model (iMantovani et al\ . \2004 ). The thin line is 
the result of a crust model with uniform density spread over a spherical shell of thickness 30 km with the detector at its vertex, 
Eq. (|B2)l . 
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FIG. 36 Uranium antin eutrino flux from the mantle as function of the distance from the detector. Points corresponds the 
25-shell reference model (iDziewonski and Andersonl . Il98ll ). The thin line describes a uniform uranium density model with the 
same mass, with uniform uranium density pu = 5.16 x 10~* g cm~'^ from a depth of 32 km up to one of 3478 km: this can 
be obtained with an average density p = 5 g cm~^ and an average abundance a\j = 10.32 ppb. The thick gray line shows a 
two-shell model, with uranium density 2.29 x 10~* g cm~^ in the upper mantle and 6.48 x 10~* g cm~"^ in the lower mantle: this 
can be obtained with average densities 3.5 g cm~^ and 5 g cm~^ and average abundances au ~ 6.55 ppb and au = 12.96 ppb, 
respectively. 
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not only on the measurements of the individual samples, but also on how the samples are chosen to be representative 
of a much larger portion of material. 

A comprehensive critical compilation of data and estimation of errors is needed. 

We remark that, when discussing uncertainties of crust and upper mantle abundances, in the past our group followed 
a different procedure that should be upgraded: we have considered the spread of published data as a full-range error, 
equivalent in same sense to a 3cr error. The numerical difference between the two procedures can be appreciated by 
looking at Table IIXI Given the small number of independent published data, this procedure underestimates the error 
compared to the standard deviation divided by \/n — 1: the case of only two data is paradigmatic. Moreover, the 
errors can only become larger the more data become available. 



2. Global and local source distributions: errors on theoretical hypotheses. 

In several instances, sought-after quantities depend on theoretical hypotheses or unknown parameters, e.g., distribu- 
tion of elements in the mantle (see Sect. lIV.Dj) or regional sources as a subducting slab under Japan (see Sect. |VlLB]). 
In such cases, we suggest to consider the minimal and maximal result as extremes of the ± 3a interval, so that the 
± 3f7 interval cover the complete range of theoretical hypotheses. 

Our group has basically followed the above-described procedure using full-range errors that included the spread of 
theoretical hypotheses. 



3. Combining errors: correlations. 

When building a reference model for geo-neutrinos, there are several sources of uncertainties (abundances, source 
distributions, oscillation parameters, cross sections, etc.) and one needs to combine them to obtain the total error. 
When errors are independent, one can combine them quadratically, but when they are correlated (anti-correlated), 
this procedure underestimates (overestimates) the total error. Correlation between U and Th abundances constitutes 
an important example: the ratio of these abundances has a smaller fluctuation than the individual abundances, 
indicating strong correlation. In fact, the only other important (anti)correlation is the one between abundances in 
the crust and in the mantle, if we impose that the total element mass (in the mantle and crust) be constrained to 
either the BSE value or to a given value used as free parameter. 

In principle it would be nice to have not only good determination of errors on the abundances, but also of their 
correlations. Given the present situation, we believe that a correct and robust approach is to select data and use them 
in such a way as to reduce correlation. 

As an example the U/Th correlation can be tackled by performing the complete calculation for uranium neutrinos 
and than scaling the result for U+ Th using a fixed ratio of U/Th: this procedure is equivalent to a 100% correlation 
(correlation coefficient p = 1) and sl ightly overestima t es the error. 

This approach was used in Ref. (jMantovani et "oZI. I2004D and also in this paper. In fact if one uses p — 0.94 as 
in Ref. (' Fogli et all l2005l ) the part of the error on the total signal due to crust abundances changes by about 1%. 
For instance the total error of our theoretical prediction at Kamioka is 5.9 TNU: if this error were due only to crust 
abundances, the use oi p — 0.94 instead of p = 1 would reduce the error to 5.85 TNU: the error itself is not known 
with that precision. 

The correlation between mantle and crust abundances is avoided not using mantle abundances as variables, but 
instead total (crust plus mantle) masses. 



